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CHAPTER  I 


INTRODUCTION 
Scope  and  Purpose 

This  dissertation  systematically  develops  estimates 
of  demand  functions  for  barge  transportation  for  a  number 
of  commodities.  The  barge  demand  functions  reflect  compe¬ 
tition  from  a  full  range  of  alternative  transport  modes. 

At  issue  are  ways  of  collecting  data  for  the  full  range 
of  attributes  which  have  a  significant  impact  on  the  choic 
of  transport  mode,  the  analysis  of  the  attributes  of  the 
data,  the  formulation  of  a  theoretically  sufficient  model 
of  transport  mode  selection  by  transport  shippers  and  im¬ 
plementation  of  the  empirical  analysis  to  generate  demand 
functions.  The  analysis  is  framed  in  the  short  run  be¬ 
havior  of  transport  users  and  thus  omits  response  which 
includes  the  relocation  of  economic  activity. 

The  analysis  is  performed  for  two  large  regions  of 
the  nation.  Because  of  the  large  amount  of  data  represent: 
ing  several  transport  modes  and  commodities,  determination 
of  the  proper  level  of  disaggregation  can  be  evaluated. 

The  model  developed  stems  from  an  extensive  investigation 
of  models  using  disaggregate  modal  choice.  The  analysis 
can  replicate  individual  shipper  mode  selection  decision¬ 
making,  taking  into  account  variables,  identifiable  in 


r 


cost  accounts  or  not,  which  influence  the  profitability  of 
a  business  firm.  If  these  variables  are  not  represented 
by  economic  cost  estimates  relevant  to  each  shipment,  the 
model  accepts  non-price  variables  and  imputes  an  economic 
weight  to  those  variables. 

Policy  Issues  and  Background 

Demand  estimates  of  barge  transportation  have  become 
more  significant  to  policy  and  project  analysis  as  Federal 
policy  interests  have  become  sensitized  to  intermodal  com¬ 
petition  for  traffic  and  Federal  subsidies.  Federal  policy 
interest  in  intermodal  competition  among  freight  haulers 
has  increased  sharply  in  response  to  the  decreasing  profit¬ 
ability  of  the  rail  mode.  Bankruptcy  proceedings,  proposals 
to  close  unprofitable  rail  line  segments,  requests  for 
massive  subsidies,  and  many  other  problems  have  become 
prominent  influences  on  federal  policy  toward  railroads. 
Important  for  the  decreasing  rail  profitability  explanation 
is  the  claim  that  federal  barge  transportation  investment 
policy  results  in  subsidies  to  the  barge  mode  that  allows 
it  to  attract  traffic  and  potential  profit  away  from  the 
competing  rail  mode. 

The  U.S.  Army  Corps  of  Engineers'  proposal  to  replace 


Lock  and  Dam  26  on  the  Mississippi  River  near  Alton,  Illi¬ 
nois,  is  a  subject  of  a  suit  brought  to  court  by  Western 


Railroads  and  environmental  interests  in  1972.  A  Federal 
court  ruled  that  authority  for  lock  and  dam  replacement 
of  the  extent  proposed  by  the  Corps,  was  insufficient  and 
thereby  restrained  construction  until  Congress  explicitly 
authorized  the  project.^  One  prominent  point  in  the  suit 
was  the  claim  that  additional  barge  capacity  was  not  needed 
because  of  excess  capacity  in  competing  rail  lines  and 
that  extra  barge  mode  system  capacity  would  adversely 
affect  rail  profitability.  One  of  the  political  reper¬ 
cussions  of  this  claim  was  a  proposal  for  a  barge  user 
tax  coupled  with  the  appropriations  action  for  Lock  and 
Dam  26 . ^ 


^The  Rivers  &  Harbors  Act  of  1909  authorizes  the 
replacement  of  existing  projects  as  a  major  rehabilitation 
action.  The  Corps  replaced  38  low  head  navigational  locks 
and  dams  on  the  Ohio  River  with  18  high  lift  structures 
and  changed  the  lock  size  from  110'  x  600'  to  twin  110'  x 
1200'  and  110'  x  600'.  These  changes  increased  the  capac¬ 
ity  of  the  system  and  reduced  the  waiting  time  for  tows 
and  maintenance  costs.  However,  the  U.S.  District  Court 
for  the  District  of  Columbia  (Civil  Action  741190,  Atchi¬ 
son,  Topeka  and  Santa  Fe  et  al  vs.  Howard  Calloway  et  al 
and  Civil  Action  741191,  Isaac  Walton  League  of  America 
et  al  vs.  Howard  Calloway  et  al)  ruled  in  the  Lock  and 
Dam  26  case  that  the  Corps  had  exceeded  its  replacement 
authority  and  that  explicit  Congressional  authorization 
for  a  new  project  was  required. 

2 

An  extensive  series  of  newspaper  articles  have  ac¬ 
companied  the  Lock  and  Dam  26  case  and  legislation  requir¬ 
ing  waterway  user  fees,  R.T.  Reid,  "A  Deft  Maneuver  by 
Barge  Fee  Backers,"  The  Washington  Post,  Tuesday,  May  2, 
1978,  p.  A4 .  Mr.  Reid  wrote  3o  articles  on  this  issue  be¬ 
tween  March  30,  1977  and  October  14,  1978. 
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If  the  topology  of  the  demand  function  for  water 
transportation  is  known,  the  benefits  of  adding  new  water¬ 
way  capacity,  the  impacts  of  imposition  of  new  policies 
for  cost  sharing,  and  the  impacts  of  changes  in  operations 
policies  can  be  determined.  Figure  1.1  shows  the  amount 
that  waterway  users  would  be  willing  to  pay  for  increased 
capacity  that  permits  additional  traffic  to  move  at  lower 
costs . 

Barge  Price 


Quantity  Shipped  by  Barge 

Figure  1.1  A  shift  in  the  supply  function  for 
barge  transportation  resulting  from 
improvements  in  operating  conditions. 

An  increase  in  the  supply  schedule  S  to  S^,  results 
in  a  decline  in  price  from  P  to  P^  and  an  increase  in 
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quantity  from  Q  to  Q^.  This  model  shows  a  simple  case 
in  which  the  added  capacity  does  not  induce  a  shift  in 
the  demand  schedule.  This  supply  schedule  is  of  the  tow¬ 
ing  industry  and  does  not  reflect  the  costs  of  altering 
the  waterway  system  to  provide  additional  capacity.  To 
accomplish  a  benefit  cost  analysis  the  area  Q  AB ,  which 
reflects  the  amount  that  transport  users  would  be  willing 
to  pay  for  the  added  capacity,  is  compared  to  the  costs 
of  adding  the  capacity.  The  difference  between  and  P 
multiplied  by  quantity  Q  is  cost  savings  to  traffic  using 
the  waterway  in  the  "without  project"  condition.  Area 
A  C  B  is  a  consumer  surplus  for  users  which  ship  addition¬ 
al  quantities  as  a  result  of  lowered  barge  transport  costs 

3 

in  the  "with  project"  condition. 


Current  Benefit  Evaluation  Requirements 
The  rules  for  federal  evaluation  of  a  navigation 

4 

project  benefits  are  stated  in  the  Procedures  Manual. 


^Charles  W.  Howe,  et  al..  Inland  Waterway  Transporta¬ 
tion,  Studies  in  Public  &  Private  Management  and  Investment 
Decisions ,  Johns  Hopkins  Press  (1969),  presents  a  much  more 
extensive  discussion  of  the  economic  issues  embodied  in  the 
analysis  of  public  and  private  investment  in  shallow-water 
inland  transportation.  Effective  public  decision-making 
regarding  investment  and  operating  policies  must  take  into 
account  the  response  of  private  carriers  and  similarly  pri¬ 
vate  carrier  decisions  must  take  into  account  present  and 
anticipated  future  public  policies. 

^U.S.  Water  Resources  Council,  "Procedures  for  Evalua¬ 
tion  of  National  Economic  Development  (NED)  Benefits  and 
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The  manual  defines  four  categories  of  benefits  of  a  naviga¬ 
tion  project: 


a.  Cost  Reduction  Benefit  (Same  Origin-Destination; 
Same  Mode) .  The  same  amount  of  traffic  would  use  the  water¬ 
way  both  "with"  and  "without  project.”  The  measure  of 
benefit  is  the  reduction  in  costs  on  an  existing  waterway 
where  improved  operating  conditions  occur  with  the  project. 


PRICE  D 

BARGE  barge 


Figure  1.2  Cost  Reduction  Benefit  Case 

A  change  in  costs  from  to  would  result  in  bene¬ 
fits  represented  by  the  hatched  area  in  Figure  1.2.  The 
following  is  a  list  of  cost  reduction  possibilities: 

(1)  Reduction  in  total  trip  delays 

(2)  Increases  in  efficiency  due  to  improvements 
which  allow  longer  or  larger  tows,  and  result  in  decreases 
in  transportation  costs  per  unit  of  goods  hauled. 

Costs  in  Water  Resources  Planning  (Level  C) ;  Proposed  Rules 
and  Procedures  of  Implementation,"  Federal  Register,  44; 

No.  102  (Thursday,  May  24,  1979),  see  paragraph  704.127, 
pp.  30220-30225. 


(3)  Increased  efficiency  due  to  barges  being 
more  fully  loaded,  as  a  result  of  deeper  channel  depths. 

b.  Shift  of  Mode  Benefit  (Same  Origin-Destination, 
Different  Mode) .  In  this  case  shippers  use  the  waterway 
"with"  the  project  but  another  mode  "without"  the  project. 
The  problem  is  to  estimate  diversion  from  one  mode  to 
another.  Impacts,  if  any,  on  barge  costs  would  derive 
from  the  impacts  discussed  in  the  cost  reduction  benefit 
case  above.  Substitutability  between  modes  is  required 
for  existence  of  this  set  of  benefits. 

c.  Shift  of  Origin-Destination  Benefit  (Different 
Mode,  Different  Origin-Destination).  In  this  case  the 
location  of  the  origin  or  destination  and  the  mode  of  trans¬ 
portation  is  different  "with"  the  project  than  "without" 
the  project.  Reductions  in  costs  other  than  transportation 
costs  may  be  included  in  the  transportation  benefit.  The 
analytical  requirements  for  this  case  are  much  more  exten¬ 
sive  than  for  the  cases  where  origin  and  destination  do 

not  change  as  a  result  of  the  project.  Because  lowered 
transport  costs  could  shift  comparative  advantage  to  an 
area  with  lower  costs  for  factors  other  than  transportation, 
benefits  estimated  may  include  the  new  location  advantages 
for  all  factors. 

d.  New  Movement  Benefit.  This  case  results  from 
new  traffic  due  to  reductions  in  transport  cost  which 
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permit  a  commodity  to  be  produced  and  transported  in  the 
"with"  project  and  not  in  the  "without"  project  case.  In¬ 
creases  in  production  and  consumption  provide  the  basis 
for  new  movements,  and  the  benefit  is  defined  to  include 
increases  in  producer  and  consumer  surplus. 

For  each  of  these  benefit  categories,  a  special 
assumption  about  the  nature  of  the  transportation  demand 
function  is  made: 

a.  Cost  Reduction  Benefit.  It  is  assumed  that  an 
absolutely  price  inelastic  demand  function  for  a  given 
transport  mode  exists. 

b.  Shift  of  Mode  Benefit.  It  is  assumed  that  a 
demand  function  exists  which  has  price  elastic  properties. 
Substitution  between  modes  is  the  key  analytical  element 
in  estimating  the  demand  function. 

c.  Shift  of  Origin-Destination  Benefit.  It  is 
assumed  that  a  price  elastic  demand  function  exists  with 
location  economies  the  key  element  in  the  benefit  set. 

d.  New  Movement  Benefit.  It  is  assumed  that  the 
transportation  demand  function  may  shift  to  the  right  due 
to  increased  goods  consumption  (with  induced  increases 

in  transportation  demand) ,  or  to  increased  production  re¬ 
sulting  from  decreased  transportation  costs  and  subsequent 
use  of  the  barge  mode. 
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T1  i  Manual  defines  Che  steps  required  for  benefit 
estimate.  "Without"  project  conditions  are  defined  to 
be  those  most  likely  conditions  expected  to  exist  in  the 
future  in  the  absence  of  a  new  project  and  any  associated 
changes  in  public  policy  towards  barge  transportation, 
including : 

(1)  All  economically  feasible  moves  within  the 
discretion  of  the  operating  agency  will  be  implemented 

at  the  appropriate  time  (including  helper  boats  and  lock 
operating  policies  which  reduce  delays) . 

(2)  User  charges  provided  by  law  at  the  time 
of  estimation. 

(3)  Normal  operation  and  maintenance  work  will 
be  performed. 

(4)  Alternative  modes  will  have  infinite  capac¬ 
ity  at  current  costs. 

(5)  Only  waterway  investments  currently  in  place 
or  under  construction  are  assumed  to  be  in  place  over  the 
period  of  estimation. 

"With"  project  conditions  are  defined  as  the  most 
likely  conditions  expected  to  exist  in  the  future  with 
the  new  project  in  place,  including: 


(1)  Demand  management  measures,  including  the 


use  of  congestion  or  lockage  fees,  that  are  to  be  evaluated 
as  nonstructural  alternatives.^ 

(2)  Actions  which  increase  the  supply  capacity 
such  as  traffic  management  are  to  be  evaluated  as  non¬ 
structural  alternatives. 

(3)  Physical  changes  which  increase  supply 
capacity,  such  as  a  new  waterway,  larger  locks,  increasing 
channel  depth,  improving  channel  geometry  in  curves  are 

to  be  evaluated  as  structural  alternatives. 

(4)  Combinations  of  nonstructural  and  structural 
alternatives . 

(5)  Alternative  timing  of  structural  measures. 

(6)  Inclusion  of  authorized  additions  in  the 
waterway  system. 

Procedures  to  be  followed  in  navigation  benefit  esti¬ 
mation  incorporate  ten  steps: 


(1) 

Identify 

commodity  types. 

(2) 

Identify 

study  area. 

(3) 

De  termine 

current  commodity  flows 

^The  manual  classifies  alternative  actions  which 
reduce  the  cost  of  transportation  as  structural  or  non¬ 
structural.  Structural  measures  include  traditional  means 
of  increasing  the  physical  capacity  of  a  waterway,  such 
as  larger  locks,  deeper  or  straighter  channels.  Non¬ 
structural  measures  tend  to  be  nontraditional  means  for 
decreasing  costs,  which  might  include  the  imposition  of 
congestion  fees  or  provision  of  helper  boats  to  reduce 
the  waiting  times  at  congested  locks. 


(4)  Determine  current  cost  of  moving  commodity 


flows  by  competing  modes. 

(5)  Forecast  commoditiy  demands. 

(6)  Project  cost  of  competing  modes. 

(7)  Determine  current  cost  of  waterway  use. 

(8)  Determine  future  cost  of  waterway  use. 

(9)  Determine  waterway  use  with  and  without 

proj  ect . 

(10)  Compute  annual  benefits  and  discount. 

Steps  9  and  10  require  an  approximation  of  or 
assumption  about  the  transportation  demand  function  for 
each  commodity  and  each  origin-destination  pair.  In 
addition,  the  choice  of  procedures  used  to  estimate  trans¬ 
port  demand  can  influence  the  estimate  of  benefits  substan¬ 
tially.  Indeed,  estimation  procedures  which  omit  transport 
demand  determinants  and  favor  supply  determinants  can  cause 
substantial  errors  in  estimating  benefits  in  every  case 
except  the  cost  reduction  benefit.  The  historic  emphasis 
on  supply-oriented  benefit  methodology  comes  from  the  fact 
that  estimates  of  product  and  transportation  demand  sched¬ 
ule  are  very  difficult  to  measure  empirically,  even  with 
good  data  representing  demand  determinants.  Coupled  with 
the  difficulties  of  estimating  transport  demand  is  the 
fact  that  transport  demand  data  are  difficult  for  public 
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agencies  to  obtain  from  privately  owned  shippers  or  trans¬ 
port  companies.^ 

In  this  dissertation  the  demand  for  waterway  trans¬ 
portation  is  estimated  from  data  obtained  from  transport 
users  selecting  between  competing  modes.  The  data  are 
used  to  form  a  mode  choice  model  which  replicates  transport 
user  behavior.  Profit  maximizing  behavior  is  assumed 
and  nonpriced  elements  of  transportation  choice  are 
allowed  to  affect  users  internal  economies.  Data  elements 
include  the  time  of  transit,  delivery  time,  handling  costs 

^The  U.S.  Bureau  of  the  Budget  issued  evaluation 
criteria  in  the  early  1960s  which  required  benefits  of 
a  waterway  project  to  be  derived  from  estimates  of  savings 
in  the  carrier  costs  rather  than  savings  to  shippers  who 
would  use  a  waterway  in  the  "with  project"  but  another 
mode  in  the  "without  project"  condition.  This  led  to 
research  conducted  by  Northwestern  University  for  the 
U.S.  Army  Corps  of  Engineers,  Leon  Moses  and  Lester  Lave 
(Eds.),  Cost  3enefit  Analysis  for  Inland  Navigation  3  Vols., 
IWR  Report  70-1  (U.S.  Army  Engineer  Institute  for  Water  Re¬ 
sources,  Ft.  Belvoir,  Va . ,  1970).  The  report  developed 
an  engineering  process  cost  function  for  barge  and  rail 
modes.  The  cost  functions  gave  a  good  approximation  of 
line  haul  costs.  However,  indirect  costs  are  not  amenable 
to  estimation  by  this  procedure  and  in  the  case  of  rail 
and  pipeline  modes,  decreasing  marginal  costs  are  held 
to  be  typical  at  least  until  capacity  limits  are  ap¬ 
proached.  These  indirect  costs  are  significant  increments 
of  total  costs.  The  report  evaluated  the  use  of  statisti¬ 
cal  cost  functions  to  estimate  indirect  costs  and  recom¬ 
mended  further  development  of  an  integrated  Engineering- 
Statistical  approach.  Moses  and  Lave  advanced  the  modal 
split  model  based  on  discriminant  analysis  and  oriented 
to  estimating  the  demand  function  for  barge  transporta¬ 
tion  as  a  superior  alternative.  However,  by  the  time 
that  the  research  was  underway,  the  U.S.  Congress  inter¬ 
vened  with  a  definition  of  benefits  in  the  1966  Transporta¬ 
tion  Act  (PL  89-670).  That  definition  holds  today,  namely 
that  benefits  should  be  measured  by  savings  to  transport 
users  based  on  rates  at  the  time  of  the  study. 


encumbered,  and  per  mile  transportation  cost.  Basic 
attributes  of  commodity  shipments  such  as  the  amount 
shipped  each  year  and  the  amount  of  each  shipment  are 
included.  Using  these  data  from  actual  shipments  and 
records,  rather  than  asking  a  series  of  "what  if  questions 
of  the  shippers,  an  estimate  of  transport  demand  is  gener¬ 
ated  which  is  free  of  the  bias  expectable  when  potential 
waterway  users  are  asked  whether  a  new  (or  improved)  water¬ 
way  would  enhance  their  profits  and  they  would  overstate 
potential  traffic  and  savings  which  could  accrue  from 
a  potential  project. 

A  general  case  is  presented  here,  for  the  total 
transportation  system.  Assume  two  transport  modes,  rail 
and  barge,  and  assume  that  the  quantity  demanded  to  be 
shipped  by  each  mode  is  a  function  of  the  relative  price 
and  service  attributes  of  each  mode: 
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Traffic  shipped  by  rail  or  barge  transport 
Price  of  rail  or  barge  transport 
Service  attributes  of  rail  or  barge  trans¬ 
port 

Technology  factors  (size  of  shipment, 
annual  shipments) 


The  demand  for  each  mode  is  limited  by  the  tech¬ 


nology  which  limits  mode  substitutability.  For  example, 
a  fully  loaded  barge  hauls  1,000  to  1,500  tons  and  most 
frequently  hauls  one  commodity  in  each  shipment.  Shippers 
of  40-100  ton  quantities  per  shipment  would  find  rail  the 
only  feasible  mode.  Thus,  the  demand  function  for  the 
rail  mode  would  have  a  segment  relatively  price  inelastic 
representing  small  shippers  demand  and  the  demand  func¬ 
tion  would  become  more  price  elastic  as  shipment  size 
approaches  barge  load  size.  When  an  expansion  in  waterway 
capacity  reduces  barge  transport  rates,  goods  previously 
shipped  by  rail  would  be  diverted  to  barges.  The  size 
of  the  diversion  would  depend  on  the  importance  of  relative 
service  attributes  of  each  mode,  and  the  ability  of  former 
rail  users  to  handle  larger  shipment  sizes. 

Organization  of  the  Study 

In  Chapter  Two,  the  theory  of  derived  demand  is 
presented  as  in  Marshall. ^  This  is  followed  with  the 
theory  of  derived  demand  for  transportation  developed 

Q 

by  Moses  and  Lave,  and  extended  to  include  influence  of 
time  in  transit,  perishability,  and  uncertainty  in  market 

^Alfred  Marshall,  Principles  of  Economics,  9th 
Edition  (McMillan  and  Company,  Ltd.,  New  York,  1901), 
pp.  381-384. 
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prices  developed  by  Allen,  and  inventory  holding  developed 
by  Baumol . ^ 

Chapter  Three  reviews  the  competing  procedures  for 
estimating  transportation  demand,  including  the  gravity 
models  in  several  specifications,  the  abstract  mode  model, 
regression  models  utilizing  aggregate  quantity  and  price 
observations  across  regions  and  commodity  groups,  and 
a  number  of  disaggregate  mode  choice  models. 

Chapter  Four  discusses  the  econometric  assumptions 
and  problems  contained  in  each  of  the  basic  procedures. 
Three  generic  procedures  are  used  in  calculating  the 
various  demand  models-,  linear  regression,  nonlinear  re¬ 
gression,  and  linear  discriminant  analysis.  Linear  dis¬ 
criminant  analysis  of  the  two  group  case  is  shown  to  be 
a  special  case  of  linear  regression. 

Chapter  Five  presents  the  estimates  of  water  trans¬ 
portation  demand  and  analysis  using  linear  discriminant 
analysis.  Transformation  of  the  data  to  a  log  normal 
distribution  provided  a  means  for  developing  the  modal 
split  model  and  ultimately  the  simulation  of  the  demand 
surface.  The  data  were  disaggregated  to  major  destination 

g 

William  Bruce  Allen,  A  Model  for  the  Demand  of 
Transportation;  The  Case  of  Air  Freight  (Northwestern 
University,  Evanston,  Illinois,  1979),  Dissertation, 
pp.  18-35. 

■^W.  J.  Baumol,  et  al,  Studies  on  the  Demand  of 
Freight  Transportation  (Mathematica ,  Inc.,  Princeton  ,N  .  J  . , 
1967),  Vol .  1. 


and  major  commodity  groups.  This  disaggregation  improved 
the  effectiveness  of  the  estimating  model. 


Chapter  Six  presents  the  summary  and  conclusions 
drawn  from  the  analysis.  The  economic  interpretation 
of  the  demand  model  and  the  potential  importance  of  shipper 
behavioral  decision-making  in  transport  demand  estimates 
are  presented. 

Appendix  A  describes  and  criticizes  the  data 
utilized  in  the  empirical  estimates  of  demand  for  water 
transportation.  The  data  utilized  represent  815  indivi¬ 
dual  shipments  totaling  122  million  tons  per  year,  with 
14  major  commodity  groups  shipped  by  nine  transport  modes. 
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CHAPTER  II 


THE  THEORY  OF  DEMAND  FOR  FREIGHT  TRANSPORTATION 

Introduction 

This  chapter  reviews  the  theory  of  derived  demand 
developed  by  Marshall,  then  the  theory  of  derived  demand 
for  transportation  developed  by  Moses,  Lave  and  Allen, 
and  finally,  the  demand  for  transportation  in  the  context 
of  an  inventory  theoretic  cost  function  as  developed  by 
Baumol  and  others.  Choice  of  transport  mode  is  accommo¬ 
dated  by  the  model.  The  simplest  case  is  based  on  trans¬ 
port  rates,  then  additional  transport  choice  variables 
are  added  (delivery  time  and  perishability) .  Next  un¬ 
certainty  in  market  prices,  delivery  time  and  damage  due 
to  perishability  are  added.  The  inventory  theoretic  model 
introduces  speculative  inventory  holding  due  to  uncertainty 
in  market  price,  delivery  time  and  perishability.  The 
discussion  is  limited  to  behavior  of  the  profit  maximizing 
firm  in  the  short  run. 

Extensive  literature  in  transport  planning  presents 
the  choice  of  transportation  problem  in  the  context  of 
physical  flows.  Economic  variables  are  noticeably  absent 
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in  most  gravity  models.^  However,  this  dissertation 
emphasizes  transport  demand  to  reflect  economic  choices 
by  managers  who  decide  how  much  to  produce  and  to  ship 
to  any  destination. 


Derived  Demand 


Transportation,  like  other  inputs  to  production, 
is  characterized  by  demand  which  is  indirect  and  derived 
from  demand  for  those  products  which  it  helps  to  produce. 
Following  is  the  classic  definition  of  derived  demand 
by  Marshall: 

The  price  that  will  be  offered  for  anything  used  in 
producing  a  commodity,  for  each  separate  unit  of  the 
commodity,  is  limited  by  the  excess  of  the  price  at 
which  that  amount  can  find  purchasers,  over  the  sum 
of  the  prices  at  which  the  corresponding  supplies 
of  the  other  things  needed  for  making  it  will  be 
forthcoming . 

To  use  technical  terms ,  the  demand  schedule  for  any 
factor  of  production  of  a  commodity  can  be  derived 
from  that  for  the  commodity  by  subtracting  from  the 
demand  price  of  each  separate  amount  of  the  commodity 
the  supply  prices  for  corresponding  amounts  of  the 
other  factors. 2 


John  R.  Meyer,  and  Mahlon  R.  Straszheim,  Techniques 
of  Transport  Planning,  Pricing  &  Project  Evaluation  (Brook¬ 
ings  Institute,  1971),  Vol.  1,  Chapter  9. 

Marshall,  Principles,  see  Chapter  6  of  Book  V, 

p.  383. 
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PRICE 


Figure  2.1  Demand  for  Knife  Handles  Based  on  the 
Demand  for  Knives  and  the  Supply 
Schedules  for  Knives  and  Blades 

The  problem  is  to  determine  the  demand  for  knife 
handles,  given  the  demand  schedule  for  knives  (D  D^)  and 
the  supply  schedules  for  knives  (S  S^)  and  handles  (S'S^) 
as  shown  in  Figure  2.1.  The  equilibrium  price  for  knives 
is  OA  at  quantity  OB.  The  demand  for  handles  (D'D^)  can 
be  derived  by  subtracting  the  differences  in  supply  prices 
of  knives  and  handles  from  the  demand  schedule  for  knives 
(D  Dk).  Let  M  ?l  cut  S'Sh  at  Mq  and  SSk  at  MQ.  The  supply 
price  of  handles  is  Oq  and  the  supply  price  of  knives  is 
OQ.  Subtracting  Qq  from  OP^  gives  OP  which  is  the  demand 
price  for  OM  handles,  D'D^  is  l°cus  °f  such  points. 
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Derived  Demand  for  Transportation 


Following  Moses  and  Lave,  the  demand  of  a  firm  for 

transportation  is  derived  from  the  demand  for  the  final 

products  of  the  firm  and  the  supply  schedule  for  all 

3 

production  costs  except  transportation.  Demand  and 
production  functions  are  given  and  the  firm  produces  one 
product  and  operates  in  a  perfectly  competitive  environ¬ 
ment  . 


7T  -  (P  -  t)  Q  -  f  (Q)  (2.1) 

where 

rr  =  producing  firm's  profit, 

P  =  market  price  of  the  firm's  product, 

Q  =  quantity  the  firm  produces  and  ships,  and 
t  =  transportation  charge. 

The  first  term  of  equation  2.1  represents  net  revenue, 
the  second  is  the  cost  of  production. 


u  =  (P 

-  t)  -  f'  (Q)  -  0 

(2.2) 

and  P  - 

t  =  f ' (Q) 

(2.3) 

with  f" 

(Q)  >  0. 

(2.  A) 

3 

Moses  and  Lave,  Cost  Benefit  Analysis,  see  Part  1. 
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Equation  2.3  states  the  usual  profit  maximizing  con¬ 
dition  that  output  is  determined  by  the  intersection  of 
marginal  cost  and  marginal  revenue  (P)  at  a  constant  price 
over  all  output  reflecting  a  competitive  market.  The 
equation  also  yields  the  firm's  demand  for  transportation 
with  t  variable  and  P  constant. 


Figure  2.2  shows  the  relationship  between  the  short 
run  supply  schedule  (MG)  of  the  firm  producing  one  product, 
the  market  price  (P  )  and  derived  transport  demand  schedule 


O  Qj  Q2  Q™  Q  0  °2  °m 


QUANTITY  PRODUCED  QUANTITY  SHIPPED 

Figure  2.2  Derived  Demand  for  Transportation 
by  Firm  Given  Average  and  Marginal 
Cost  Schedules  and  Market  Price 

While  the  average  variable  cost  curve  (AVC)  is  U 
shaped,  the  demand  function  is  truncated  but  has  the  usual 


The  truncation  occurs  because  the  firm 


negative  slope. ^ 
will  shut  down  production  of  P-t^  <  since  transport 
costs  would  reduce  marginal  revenue  to  less  than  average 
variable  cost.  The  diagrams  have  illustrated  the  case 
that  at  rate  t^,  is  produced  and  shipped,  at  lower 
rates  greater  quantities  are  shipped  but  at  higher  rates 
nothing  is  produced  and  shipped. 

Additional  Transportation  Demand  Variables 


The  above  formulation  can  be  extended  to  introduce 
additional  transportation  choice  variables.  Positive 
costs  for  time  of  transit  can  be  added  as  follows:^ 

tt  =  -  f (Q)  (2.5) 

(l+i)a 


where  i 
a 

3  tt 

and, 


=  interest  rate  per  unit  of  time  and, 
=  time  of  transit. 

=  —  -  f '  (Q)  =  0 

(l+i)a 


(2.6) 


-P~—  =  f '  (Q)  ,  (2.7) 

(l+i)a 

with  f " (Q)  >  0.  (2.7a) 

^ Ibid ,  p .  3  . 

^Allen,  A  Model  of  the  Demand  for  Transportation, 
Chapter  2. 
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For  equation  2.5  it  is  assumed  that  costs  of  produc¬ 
tion  are  incurred  instantly,  and  that  the  total  time  hori¬ 
zon  is  one  time  period,  therefore  the  total  product  is 
sold  instantaneously  once  the  goods  arrive  at  the  market 
and  payment  for  the  product  is  COD.  There  is  a  time 
separation  between  production  of  goods  and  the  receipt 
of  revenues  from  their  sale.  In  this  equation,  the  second 
term  is  the  total  cost  of  production  and  the  first  term 
is  the  present  value  of  the  revenue.  The  derivation  of 
the  demand  for  transportation  follows  exactly  in  the  same 
fashion  as  developed  on  page  20.  If  the  average  cost 
function  is  U  shaped,  the  transport  demand  function  is 
truncated  but  has  the  usual  negative  slope.  The  trunca¬ 
tion  occurs  because  the  firm  will  stop  production  if  the 
transport  cost  rises  to  the  point  where 

(P-t)  =  p 

(l+i)a 

is  less  than  minimum  average  variable  cost. 

Allen  introduced  a  term  for  perishability  as  follows 

tt  =  (1~B)  (P"^  -  f (Q)  (2.8) 

( 1+i) a 

Where  B  =  damage,  pilferage  losses  or  perishability 
rate,  and  the  numerator  of  the  first  term  represents  net 
revenue  after  deduction  for  perishability.  Taking  the 


partial  derivative  of  profit  with  respect  to  quantity, 
we  have : 


£  -  -  *'<q)  -  o 

(1+i) 

(2.9) 

and 

( 1  —  B ) (P-t)  _  f , (Q) 

(2.10) 

( 1+i) a 

with 

f"(Q)  -  0. 

(2.11) 

Again  by  separating  the  time  of  production  of  the 
goods  and  realization  of  revenues,  the  model  accounts  for 
losses  in  revenue  due  to  perishability  during  transit 
(B) .  Net  revenues  reflect  both  perishability  and  costs 
of  time  during  transit.  The  transportation  demand  function 
derived  from  this  specif ication  is  truncated  and  negative¬ 
ly  sloped. 

Allen  developed  the  model  to  include  uncertainty 
in  market  prices  and  damage  rate  due  to  perishability . ^ 
Given  uncertainty  of  transport  costs  (t),  and  interest 
rates  (i),  the  firm  would  assume  an  average  market  price 
(P)  and  average  damage  rate  (B) . 

n 

Define  (P)  =  £  PkXk  w^ere  lXk^  is  C^e  Pro^a- 

k=l 

bility  distribution  of  prices. 

6Ibid. ,  p.  27. 
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n 

Define  (B)  =  Y  B  Y  where  {Y  }  is  the  proba- 

1Z  z  z 

bility  distribution  of  damage  factors. 


n  -  (1-B) (P-t)  0 

( 1+i) a 

f  (Q) 

(2.12) 

3 it  (1-B)  (P-t) 

3(5  (l+i)a 

f ' (Q)  =  0 

(2.13) 

and 

(1-B)  (P-t)  _  f.r) 

(l+i)a 

(2.14) 

With 

f ” ( Q )  >  0. 

(2.15) 

This  model  gives  the  same  general  result  as  previous¬ 
ly  discussed.  The  derived  demand  for  transportation  is 
truncated  but  negatively  sloped. 

If  delay  times  are  varied  the  analysis  is  axtered. 
Following  Allen,  instead  of  assuming  an  expected  delay 
time  prevails,  the  firm  would  maximite  profits  by  producing 
the  outputs  determined  by  the  expected  net  discounted 
price  (P)  : 

p  =  Ikil  LLiI  y  (2.16) 

(l+i)am  m 

where  ym  is  the  probability  distribution  of  delays,  and 

P  >  (P-t).  (2.17) 

(1+i) a 

because  the  denominator  is  not  linear  in  its  parameter. 
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With  all  three  variables  random,  the  net  discounted 


price  would  be 


P 


zzz 

mzk 


( 1+i) am 


(2.18) 


Assuming  speculative  inventory  holding,  the  analysis 
shifts  to  a  desire  to  hold  inventory  because  of  variable 
market  prices  {P^}.  Define  (P^l  =  w^ere  {X^} 

is  the  probability  distribution  of  market  prices  {P^}  and 
q  =  inventory  costs  per  unit  per  day. 

Simplifying  the  case  by  constraining  the  inventory 
holding  period  to  one  day,  the  firm  will  engage  in  specu¬ 
lative  inventory  holding  in  its  output  decision  if, 


( l+i) 


(2.19) 


The  above  (2.19)  states  that  if  the  discounted  value 
of  production  for  inventory  is  greater  than  current  market 
price  P-p  then  firm  will  produce  for  inventory  and  shipment 
Where  P  is  the  expected  value  of  market  prices  {P^} 
and  q  is  discounted,  inventory  costs  and  profit  maximizing 
output  is  equal  to: 


(!?!)  X1  +V2  +P3X3 


Q  -  f(Q) 


(2.20) 
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^  _  P-q  v  .  p  v 

3q  ( 1+i)  *1  +  zA2 

+  P3X3  -  f ' (Q)  =  0 

(2.21) 

and  X,  +  P0X0  +  P-X,  = 

(l+i)  1  22  33 

f  ’  (Q) 

(2.22) 

with  f " (q)  >  0. 

(2.23) 

Baumol,  et  al,  extends  the  case  to  inventory  holding 
in  a  multi-time  period  analysis  based  on  an  inventory 
theoretic  cost  function.'7  Define  costs  (C)  equal  to  direct 
shipping  costs  plus  intratransit  carrying  costs  plus  inven¬ 
tory  carrying  costs. 


C  =  RT  +  vzT  +  wsT/2  +  a/s 


(2.24) 


Where  C 
R 
T 
v 

z 

a 

s 

w 


=  expected  total  variable  cost, 

=  transport  rate, 

=  transport  flow  per  year, 

=  carrying  cost  in  transit  (interest 
(plus  deterioration  plus  pilferage) , 

=  transit  time, 

=  order  cost, 

=  time  between  shipments, 

=  warehouse  carrying  costs  per  unit  per  year. 


The  shipper  controls  reorder  frequency  and  his  optimal 
decision  will  follow  the  rule: 


/ 


^Baumol  et  al. 
Transportation,  Vol. 


Studies  on  the  Demand  for  Freigt 
1,  Chapters  1-3. 
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therefore,  s  =  \2a/wT, 

and  substituting  for  s  in  2.24  and  simplifying 


C  =  RT  +  vzT  +  a/  yj  2a/wT  +  (wT/2)  V2a/wT  (2.26) 

=  RT  +  vzT  +  a  VwT/2a  +  (  V  wT / 2 )  2  >/2a/wT  (2.27) 

=  RT  +  vzT  +  wT/2  +  -yj  wT/2  VT  (2.28) 

=  RT  +  vzT  +  V2awT  (2.29) 


Uncertainty  of  delivery  time  and  demand  requirements 
on  inventory  suggests  a  safety  stock.  If  the  uncertainty 
is  characterized  by  a  Poisson  distribution,  safety  stock 
level  (S)  is: 

S  =  k  V (s  +  z)  T  (2.30) 

where  k  =  a  constant, 

and  y (s+z)T  =  standard  deviation  of  available  inventory ,  which 
leads  to  an  expanded  cost  equation, 


(2.31) 


In  this  case,  delivery  time  (z)  enters  the  inventory  figure, 
and  given  uncertain  timing  in  demands  for  the  commodity, 
slower  delivery  would  require  an  increase  in  safety  stock. 
Differentiating  (2.31)  with  respect  to  T  allows  the  margi¬ 
nal  shipping  costs  to  be  determined: 


3C 

3T 


R  +  V2  + 


(2.32) 


But  this  formulation  indicates  the  shipper  is  simply  try¬ 
ing  to  minimize  transport  costs  rather  than  maximizing 
profits,  so  a  revenue  term  should  be  added: 


MR  =  AP  +  t  ^  (2.33) 

Where  MR  =  marginal  revenue 

and  aP  =  price  differences  between  origin  and  desti¬ 
nation  (normally  equal  to  transportation  costs)  and  sub- 
d'  P 

stituting  b  cor  set  marginal  costs  equal  to  marginal 

revenue : 


=  ip+bT 
(2.34) 


This  results  in  a  nonlinear  equation  requiring  estimation 

O 

by  simultaneous  utility  search  techniques. 

Because  space  and  distance  do  influence  price  of 
transportation  (in  terms  of  price  per  unit  of  distance 
and  the  carrying  costs  associated  with  time  in  transit) , 
the  above  formulation  is  incomplete.  Distance  and  competi¬ 
tion  between  modes  are  introduced  in  the  next  chapter. 


g 

Ibid,  p.  53.  The  Gradx  computer  routine  developed 
by  R.E.  Quandt  is  one  of  the  developments  which  permits 
estimates  of  the  parameters . 
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CHAPTER  III 


TRANSPORT  DEMAND  MODELS 


Gravity  Models 


Numerous  studies  of  transport  flows  have  been  organ¬ 
ized  around  the  physical  concept  of  attraction  and  impe¬ 
dance.  ^ 


,  Ii!i 

lo  d .  . 


(3.1) 


Where  T^ 
j  cities, 


is  the  measure  of  traffic  between  i  and 


P^,  Pj  are  the  populations  of  cities  i  and  j, 
dij  rePresenCs  the  distance  between  i  and  j ,  and 
aQ  is  a  parameter  to  be  estimated  from  the  data. 


This  "gravity"  model  is  roughly  similar  to  the  gravi¬ 
tational  force  formula  in  Newtonian  mechanics.  The  volume 
of  traffic  between  cities  i  and  j  is  directly  related  to 


^See  Roger  E.  Alcaly,  "The  Demand  for  Air  Travel," 
in  Studies  in  Travel  Demand  (Princeton,  Mathematica,  1966), 
12;  and  "Aggregate  and  Gravity  Models:  Some  Empirical  Evi¬ 
dence,"  Journal  of  Regional  Science,  7  (1967),  pp .  61-73; 
also  Walter  Isard,  Methods  of  Regional  Analysis :  An  Intro¬ 
duction  to  Regional  Scienci  (Joseph  Wiley  and  Sons,  1960) , 
Chapter  II;  and  Richard  E.  Quandt,  "Some  Perspectives  of 
Gravity  Models,"  Studies  in  Travel  Demand  (Princeton, 
Mathematica,  1966)~  2^  Chapter  3. 
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the  product  of  the  two  populations,  i.e.,  the  attraction, 
and  inversely  related  to  the  distance  between  them,  i.e., 
the  impedance,  with  the  relation  between  attraction  and 
impedance  being  proportional. 

According  to  Isard: 

A  basic  attribute  of  it  (gravity  model)  used  for 
projections  is  the  lack  of  theory  to  explain  the 
values  or  functions  which  we  can  assign  to  weights 
of  exponents.  Currently,  the  justification  for 
the  gravity  model  is  simply  that,  everything  else 
being  equal,  the  interaction  between  any  two  popu¬ 
lations  can  be  expected  to  be  directly  related  to 
their  size,  and  since  distance  involves  friction, 
inconvenience,  and  cost,  such  interactions  can  be 
expected  to  be  inversely  related  to  distance. 2 

A  more  sophisticated  formulation  is: 


pflP|2 

*o  d.  . 

3-J 


(3.2) 


where  aQ ,  a^ ,  &2,  are  parameters  to  be  estimated  from 
the  data.  This  equation  is  not  a  simple  statement  of 
proportionality  of  product  of  population  and  distance 
but  indicates  an  attenuation  of  the  attraction.  The 
parameter  estimates  are  generated  by  regression  of  cross 
section  data.  According  to  Quandt,  such  regressions  often 
yield  high  correlation  coefficients  with  regression  co¬ 
efficients  significantly  different  from  zero,  but  frequent¬ 
ly  yield  bad  predictions  for  projected  travel  for  year  x^n 


2 


Isard,  Methods  of  Regional  Analysis, 


P- 


515. 


31 


from  the  cross  section  of  data  obtained  for  the  year  x. 

It  follows  that  the  bad  predictions  are  a  result  of  a 
fundamentally  erroneous  paradigm  or  the  result  of  specifi¬ 
cation  errors.  Quandt  proposed  a  greatly  expanded  formula¬ 
tion  to  correct  errors  of  specification  as  follows: 


T  .  kaU>p  al(x)  p  a2<x>  v  a3 (x) 
ij  ~  K  ric  rjt  yit 


m  a,-  (x)  a,  ( x)  a7(x) 

m.„  o  m .  o  n .  .  _  /  t .  .  _ 

it  jt  ljt  ljt 


a  (x) 
ag  (x) 


Sijt 


a9(x) 


(3.3) 


where  x  =  time  period, 

i,j  =  transport  nodes  (location  i,  j), 
y  =  income , 

m  =  city  type  classified  by  portion  of  labor 
force  in  mining  and  manufacturing, 

n  =  travel  time, 

t  =  travel  cost, 

s  =  transportation  supply  characteristics ,  and 

a,k  =  parameters. 

The  resulting  regression  equation  is  linear  in 
logarithms  which  allows  for  decreasing  returns  to  the 
independent  variables.  The  expanded  formulation  was  used 
in  the  analysis  of  intercity  passenger  demand  in  the 


p.  35. 


Quandt,  ’’Some  Perspectives  of  Gravity  Models,” 


i 


Northeast  Corridor  Study  for  the  U.S.  Department  of  Trans¬ 
portation.^  Meyer  &  Strazheim  argue  that  individual  ele¬ 
ments  of  travel  time  should  be  disaggregated,  but  data 
problems  were  deemed  to  be  unfeasible.^  Cross-section 
data  were  used  to  estimate  the  functions  with  a  result 
that  many  of  the  variables  were  highly  correlated. 

Limited  empirical  applications  of  the  gravity  models 
are  contained  in  the  literature.  Alcaly  fitted  gravity 
models  as  formulated  in  (3.3)  for  air  travel  in  seven 
cities  in  the  northwest  corridor  of  the  U.S.  and  found 
that  the  regression  coefficient  of  the  distance  variable 
to  be  insignificant . ^  While  the  equations  were  fitted 
to  predict  demand  values  on  given  air  travel  routes  in 
1950,  the  results  were  uniformly  low  except  for  two  cities. 
Alcaly  cites  Kessler  for  an  example  in  which  the  basic 
gravity  model  formulation  of  the  demand  for  transportation 
was  modified  to  include  various  socio-economic  and 


System  Analysis  and  Research  Corporation,  Demand 
for  Intercity  Passenger  Travel  in  the  Washington-Boston 
Corridor  (Cambridge,  Massachusetts,  System  Analysis  and 
Research  Corporation,  1963) ,  describes  the  development 
and  calibration  of  the  model.  Gerald  Kraft  and  Martin 
Wohl,  ’’Special  Survey  Paper,  New  Directions  for  Passenger 
Demand  Analysis  and  Forecasting,"  Transportation  Research, 
1  (Nov.  1967)  review  the  arguments. 

^Meyer  and  Strazheim,  Techniques  of  Transport  Plan¬ 
ning,  pp.  139-140. 

^Alcaly,  "The  Demand  for  Air  Travel,"  p.  69. 


demographic  characteristics  of  the  cities  involved,  although 
Kessler's  analysis  did  not  include  all  suggested  charac¬ 
teristics  .  ^ 

Alcaly  later  explored  the  effect  of  aggregation  on 
gravity  models  specifically  to  determine  "the  effects 
of  aggregation  over  modes  of  travel  on  the  performance  of 

Q 

equation  (3-3)  in  explaining  the  demand  for  travel."  The 
study  used  16  pairs  of  cities  in  California  and  1960  data. 
Travel  time  as  well  as  cost  were  used  in  some  of  the  re¬ 
gressions.  Taking  account  of  the  explanatory  power  and 
significance  of  the  postulated  log  linear  relationships, 
aggregated  mode  (to  total  transport)  equations  performed 
better  as  a  group  than  did  individual  mode  equations.  An 
exception  was  the  automobile  mode.  Air  mode  equations 
were  characterized  by  high  explanatory  power  and  signifi¬ 
cance  but  distance  coefficients  were  not  significant. 

The  study  concluded  that  the  gravity  model  explains  travel 
by  all  modes  better  than  travel  by  individual  modes.  It 
follows  that  the  theoretical  objections  are  minimized 
when  the  model  uses  aggregated  data. 

^D.S.  Kessler,  Relationship  Between  Inner  City  Air 
Passenger  and  Demographic  Factors  -  A  Multiple  Region 
Analysis ,  Princeton  University,  Princeton,  N.J.,  1965), 
Thesis . 

8 

Alcaly,  "Aggregation  and  Gravity  Models;  Some  Em¬ 
pirical  Evidence,"  pp .  61-73. 
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Intervening  Opportunities  Model 


The  intervening  opportunities  model  has  the  objective 

to  minimize  the  travel  time  between  every  region,  subject 

to  the  constraint  that  every  potential  destination  (inter- 

9 

vening  and  ultimate)  is  considered.  Thus  a  trip  origi¬ 
nating  in  region  i  has  less  probability  of  ending  in  region 
j  as  the  number  of  intervening  opportunities  increases. 


Tu  -  °i  D<vi>  - 


p(Tj) 


total  probability  that  a  trip  will  terminate 
before  the  j th  possible  destination  is 
considered , 


(T.)  =  possible  destinations  already  considered 

-*  that  are  reached  before  reaching  zone  j,  and 

0.  =  the  constant  probability  of  a  possible 

1  destination  being  accepted  if  considered. 

Thus,  the  number  of  trips  taken  between  i  and  j  is 


the  number  of  regions  multiplied  by  the  probability  of  the 


trip  terminating  in  j ,  or 


(3.5) 


where  L  shapes  the  distribution  of  trips  between  i  and  j . 

As  Isard  proves,  the  results  are  identical  with  the  gravity 


9 


Meyer  &  Strazheim,  Techniques 


of  Transportation 


Planning,  p.  121. 


model  when  the  number  of  intervening  opportunities  are  a 
linear  function  of  distance.^ 


Fratar  Method 


The  Fratar  method  is  a  procedure  for  averaging 
differential  growth  factors  for  existing  traffic  counts  and 
represents  a  logical  extension  from  the  simple  growth  fac¬ 
tor  model. ^  Future  interzonal  travel  forecasts  are  de¬ 
rived  from  the  present  level  of  interzonal  trips  and  the 
different  zonal  growth  factors.  Each  zone  growth  factor, 
factor  g^  is  a  ratio  of  future  to  present  trip  generation: 

8i  =  (3.6) 

1  T, 


Where  =  future  trip  generation,  and 


T.^  =  present  trip  generation. 


=  T.  . 


zonal  interchange,  T^j  , 

is  given  by 

N 

z  i  Fr : . 
j=i  l  l.i 

*  6i 

(3.7) 

N  r 

I  T.  - 

j-1  LlJ 

'  gij 

Walter  Isard,  Location  and  Space  Economy:  A 
General  Theory  Relating  to  Industrial  Location,  Market 
Areas,  Land  Use,  Travel  and  Urban  Studies  (Cambridge, MIT , 
1956) ,  pp .  540-544. 

^See  Thomas  J.  Fratar,  "Forecasting  Distribution  of 
Interzonal  Vehicular  Trips  by  Successive  Approximation," 
Proceedings  of  the  33rd  Meeting  of  the  Highway  Research 
Board  (Washington,  D.C.,  1954);  33. 
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tut 


The  model  is  used  in  a  successive  approximation  technique 
iteratively  to  update  growth  factors  based  on  trip  desti¬ 
nations  . 

All  of  the  procedures  described  above  are  dominated 

by  empirical  generalization  and  lack  a  fully  articulated 

behavioral  b^sis  with  the  exception  of  Quandt's  sophisti- 

1 2 

cated  specifications  given  in  equation  3.3  above.  They 
are  not  useful  to  derive  fully  and  implement  a  model  of 
transport  demand. 

Abstract  Mode  Model 

The  demand  models  discussed  above  define  each  mode 
by  administrative  ownership  or  type  of  physical  equipment 
utilized  and  a  single  demand  equation  is  estimated  for 
each  mode.  In  contrast,  the  abstract  mode  approach  defines 
choice  of  mode  in  one  equation.  This  approach  results 
in  significant  data  savings,  since  the  number  of  observa¬ 
tions  available  for  estimating  one  abstract  mode  equation 
will  be  equal  to  the  number  of  modes  times  the  number 
of  city  pairs  in  the  transport  system.  By  comparison 
any  formulation  requiring  an  equation  for  each  model  re¬ 
sults  in  a  sample  size  equal  to  the  number  of  city  pairs. 

Quandt  and  Baumol  formulated  the  model  of  abstract 
mode  to  model  passenger  travel  in  the  United  States 

1 2 

Quandt,  "Some  Perspectives  on  Gravity  Models," 
pp.  33-46. 
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Northeast  Corridor.  A  single  equation  represents  passen¬ 
ger  choice  for  all  modes  and  is  linear.  Each  mode  is 
characterized  by  several  independent  variables  which 
describe  its  supply  characteristics. 

The  formulation  defines  the  characteristics  for 
which  a  mode  should  be  developed  in  the  future.  With 
changing  technology  the  approach  offers  a  definite  forward 
thrust  in  analysis. 


Si  Si 

ok  ®  7  f  ^  (h)  f  9  (c)  f  3  (d) 


(3.8) 
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transport  flow  from  node  i  to  node  j  by 
mode  k , 

population  of  nodes  i  and  j  , 

median  income  of  nodes  i  and  j , 

an  index  representing  industrial  charac¬ 
teristics  of  nodes  i  and  j , 

best  available  travel  time  between  nodes 
i  and  j  , 

relative  travel  time  for  kth  mode, 

least  cost  to  travel  between  mode  i  and  j, 

relative  costs  for  kth  mode, 

best  departure  figure  between  nodes  i  and  j , 
relative  departure  frequency  for  kth  mode, 
number  of  modes  serving  nodes  i  and  j  , 
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,  and 

a,B  =  parameters- 

The  derived  model  is  linear  in  logarithms  with  attractive¬ 
ness  and  impedance  variables.  Fundamentally,  the  specifi¬ 
cation  argues  that  the  choice  of  transport  mode  is  based 
on  a  comparison  of  a  particular  mode's  performance  to 
that  of  the  best  available  mode  in  one  particular  variable 
(cost,  time,  or  other  single  service  factor).  This  leads 
to  a  discontinuity  of  preferences  with  respect  to  mode 
of  choice  since  the  improvements  of  several  characteristics 
by  any  mode,  short  of  becoming  the  best  mode  with  that 
characteristic,  have  no  effect  on  demand. 

Meyer  and  Strazheim  conclude  that  "empirical  applica¬ 
tion  of  the  abstract  mode  formulation  has  not  been  notably 
14 

successful."  The  model  was  applied  to  sixteen  city 
pairs  for  air  and  bus  and  automobile  traffic.  Quandt 
and  Baumol  concluded  that  the  model  has  promise  on  the 

^Meyer  and  Strazheim,  Techniques  of  Transport 
Planning ,  pp .  144-151. 
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basis  of  correctness  and  signs  on  regression  coefficients, 
stability  of  exponents  from  regression  and  lack  of  prob¬ 
lems  from  multicollinearity . ^  However,  Young  found  that 
the  single  mode  demand  equations,  using  the  F  test  to 
compare  residual  sums  of  squares  of  combined  regression 
and  the  sum  of  two  residual  sum  of  squares  of  the  two 

1 6 

separate  regressions  was  characterized  by  poor  results. 

One  application  was  to  reduce  "abstractness"  by 
introduction  of  dummy  variables  or  other  modifications 
in  the  functional  form.b^  The  following  equation  shows 
introduction  of  dummy  variables  for  auto  and  bus 


T...  =  a  (p.p.)  1(c..b)a2  (c,.,r)a3  (h..b)a4 

kij  rj  kij  ij 


rN  a5 


<dkiir,a6  nii“7  * 


(a 


lkij 


J-J 


ij 


ij 


8  + 


a9gkij  +  a10bkij) 


where  and  b^^j  are  dummy  variables  representing  auto 

and  bus  mode. 


^Richard  E.  Quant  and  William  J.  Baumol,  "The  Ab¬ 
stract  Mode  Model:  Theory  &  Measurement,  Studies  in 
Travel  Demand  (Mathematica ,  Princeton,  N.J.,  1967),  2, 
p"!  53"! 

bbQuandt  and  Young,  "Cross  Sectional  Travel  Demand 
Models:  Estimates  &  Tests,"  pp.  39-74. 

h^Kan  Hua  Young,  "Testing  the  Adequacy  of  the  Linear 
Abstract  Mode,"  Studies  in  Travel  Demand  (Mathematica, 
Princeton,  N.J.,  19b  / )  ,  3"!  p"!  TO. 
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Quandt  and  Young  developed  a  number  of  other  models 
to  introduce  nonabstractness.  The  hypothesis  that  is 
not  the  "best"  cost  and  journey  time  that  are  relevant 


was  developed  by  expressing  the  same  figures  relevant 
to  income  as  measured  by  c^j  and  h  k/y_  .  Other 

variations  assume  income  elasticities  along  different 


routes  to  be  different  (but  the  same  for  different  modes 


on  the  same  route) .  Another  variation  relaxes  the  assump¬ 
tion  of  equal  income  elasticities  for  different  modes 
by  introducing  dummy  variables  for  each  mode.  Further 
variations  were  added  by  assuming  that  it  is  the  ratio 
of  relative  cost,  journey  time  and  departure  frequency 
to  income  that  is  relevant  to  travel  demand.  Two  versions 


introduced  the  concepts  embodied  in  intervening  opportuni¬ 
ties  which  assume  that  an  individual  considering  a  trip 
will  consider  not  only  the  absolute  and  relative  character¬ 
istics  of  the  modes  between  i  and  j  but  also  may  include 
the  option  to  go  to  another  mode  1.  Equations  were  fitted 
using  both  the  California  data  and  data  from  a  sample' 
of  intercity  flows  in  the  Northeast  Corridor.  Significant 
differences  in  income  elasticity  among  city  pairs  were 
revealed.  Values  ranging  from  almost  .one  to  over  three 
resulted;  an  outcome  not  easy  to  interpret. 

It  may  be  concluded  that  an  abstract  mode  formulation, 
which  represents  travel  choices  among  modes  as  a  single 
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equation,  with  each  mode's  performance  reported  relative 
to  the  time  or  cost  of  the  best  mode  serving  a  city  pair, 
does  not  appear  superior  to  conventional  multiequation 
gravity  formulation.  The  use  of  the  more  disaggregate 
form  of  modal  choice  model  is  supported  by  the  arguments 
that  (1)  the  underlying  utility  function  and  the  demand 
functions  are  too  complex  to  be  represented  in  one  equa¬ 
tion  and  (2)  by  empirical  tests. 

Freight  Demand  Modeling 

The  above  models  were  developed  primarily  to  estimate 
derived  demand  for  and  flows  of  passenger  transportation , 
although  the  principles  used  have  been  extended  to  freight 
traffic  models.  Transport  flows  of  freight  used  in  pro¬ 
duction  processes  and  passengers  behave  quite  differently 
since  in  the  first  case  managers  are  maximizing  profits 
and  in  the  second  case  passengers  are  maximizing  utility, 
at  least  in  cases  other  than  business  travel.  Transporta¬ 
tion  demand  analysts  are  also  concerned  with  the  relation¬ 
ship  between  price  and  quantity  demanded  as  well  as  in 
the  maximization  of  an  objective  function. 

The  following  are  reviews  of  the  literature  of 
freight  demand  modeling  including  modal  split  modals. 
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Regression  Models 


Perle  carried  out  an  extensive  analysis  of  five 

commodity  groups  based  on  aggregate  data  from  nine  regions 

18 

of  the  United  States.  Two  demand  models  were  developed. 
The  basic  model  was: 


r,t 


=  a. 


+  aicr  +  a2cc 


(3.10) 


where  T  =  quantity  transported  by  mode  r  or  t ,  and 

r ,  c 

cr>  ct  =  price  of  transport  by  rail  (r)  or  truck  (t) 

Analysis  was  conducted  with  all  observations  aggre¬ 
gated  to  the  nation,  then  separated  by  commodity. 

Perle' s  results  show  the  expected  inverse  relation¬ 
ship  between  price  and  quantity  demanded  and  the  results 
permitted  estimates  of  the  price  elasticity  of  demand. 

His  analysis  would  permit  the  evaluation  of  regional  im¬ 
pacts  of  changes  in  rail  or  truck  rates  on  transportation 
demand  but  does  not  permit  analysis  of  the  impacts  of 
rate  changes  on  the  level  of  service,  change  in  economic 
activity  or  on  shipment  size. 

The  second  model,  which  emphasized  the  elasticity 
of  substitution  was  formulated  as: 


1 8 

Eugene  D.  Perle,  The  Demand  for  Transportation: 
Regional  and  Commodity  Studies  in  the  United  States  (The 
University  of  Chicago  Press,  Chicago,  Ill. ,  1964) . 


(3.11) 


Changes  in  the  elasticity  of  substitution  can  be  a 
result  of  the  relative  changes  in  prices  or  relative  changes 
in  quantities.  Thus,  the  first  model  related  absolute 
quantities  and  prices  and  treated  substitutability  tangen¬ 
tially.  The  second  model  focused  on  relative  quantities 
and  prices  and  therefore  treated  substitutability  directly. 
With  nine  regions,  five  commodity  groups,  and  five  years 
of  data,  Perle's  analysis  is  based  on  a  total  225  observa¬ 
tions.  The  data,  1956  through  1960,  showed  a  sharp  six 
percent  gain  in  truck  traffic  over  rail.  The  analysis 
failed  to  consider  handling  costs,  travel  time,  frequency 
of  service  or  other  factors  in  transportation  demand  and 
the  price  information  was  limited  to  published  rates, 
which  fail  to  represent  a  large  share  of  truck  movements. 

Sloss  developed  a  model  for  estimating  intercity 

freight  shipments  to  be  carried  by  truck,  using  Canadian 
19 

statistics.  His  formulation  utilized  aggregate  province 
statistics  for  intercity  truck  traffic  (T^) ,  average  truck 
revenue  per  ton  (ct)  ,  average  rail  revenue  per  ton  (cr) , 

19 

J.  Sloss,  "The  Demand  for  Intercity  Motor  Transport: 
A  Macroeconomic  Analysis,"  The  Journal  of  Business  (1971). 


44 


•tab 


and  a  variable  representing  general  economic  activity  (h) . 
The  general  economic  activity  variable  was  composed  of 
the  sum  of  farm  income,  value  of  building  permits  and 
value  of  shipments  of  manufactured  goods.  The  following 
model  was  formulated: 


1  2  1 

where  a  ,  a  and  a  are  parameters  to  be  estimated  from 
the  data. 

Sloss  found  that  the  volume  of  truck  traffic  is 
directly  related  to  the  price  of  rail  service  and  inversely 
related  to  the  price  of  truck  service.  Further,  the  quan¬ 
tity  of  truck  traffic  is  positively  correlated  with  the 
general  level  of  economic  activity.  Thus  the  model  could 
be  used  in  a  macro  level  evaluation  of  alternative  pricing 
strategies  but  does  not  include  service  attributes  such 
as  transit  time  and  frequency  of  service  or  handling  costs. 
In  addition,  it  has  limited  capacity  to  predict  the  impacts 
of  changes  in  prices  on  transport  flows  between  particular 
origins  and  destinations. 

Subsequently,  A.D.  Little  developed  an  aggregate 

level  modal  split  analysis  in  a  contract  study  for  the 

2  0 

Maritime  Administration.  The  data  was  aggregated  by  BEA 
20 

A.D.  Little,  Inc.,  "Domestic  Waterborne  Market 
Analysis  -  Appendix  C  -  Modal  Share  Estimate,"  prepared 
for  the  Maritime  Administration  (U.S.  Dept,  of  Commerce, 
Washington,  D.C.,  1974). 
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(Bureau  of  Economic  Analysis,  U.S.  Department  of  Commerce) 
economic  region  and  the  model  was  directed  toward  correla¬ 
tion  of  the  fraction  of  plants  located  on  water  bodies 
to  the  demand  for  barge  transportation.  The  model  is 
f ormulated  as : 


1  manufacturing  plants  on  water  = 


f  (T  p  b  a  iil  ill) 

VTij  ’  pc  >bij  ’  s  ij  ’  dwi j  ’  dwi j  / 


(3.13) 


where  T^ 


s  .  . 


ij 


r  ,w,y 


total  flows  between  regions  i  and  j , 

value  per  ton  of  commodity  c, 

percent  commodities  transported  in  bulk, 

seasonality  of  shipment  between  regions 
i  and  j  , 

distance  between  regions  i  and  j ,  and 
rail,  water  and  truck  modes. 


A  nonlinear  model  was  estimated  for  15  commodity 
groups  between  BEA  areas.  The  model  results  shed  consider¬ 
able  insight  into  the  interaction  of  nonrate  modal  per¬ 
formance  characteristics  with  shipper  and  market  attributes 
on  location  of  manufacturing  plants.  Since  rate  relation¬ 
ships  are  missing,  the  model  is  inadequate  to  analyze 
the  impact  of  rate  and  level  of  service  changes  on  traffic 
flows  and  is  therefore  inadequate  to  estimate  transporta¬ 
tion  demand. 
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Mathematica  (1967)  adapted  the  Quandt/Baumol  abstract 


mode  model  in  a  contract  for  the  Northwest  Corridor  pro- 
21 

ject.  The  general  form  of  the  adapted  model  is: 


where  T.  .  . 

kij 


a-i  a7  a_  a,  a  a , 
a  P,  1  P.  2  y.  3  y.  4  m.  5  m,  6 

Ol  J  J  1  J  j  L  J 

a 7  b  B,  r.B  ,  b>  B 

"u  7  >  1  ‘“ku  1  2  (ciJ  ’  3 


,  r .  B, 

<ckij  >  4 


(3.14) 


quantity  of  freight  traffic  from  i  to  j 
by  mode  k, 


P. 


population  of  origin  and  destination  loca¬ 
tions  , 


yi’Yj 


gross  regional  product  of  origin  and 
destination  regions, 


mi’mj 


indicators  of  industrial  character  of 
regions  (such  as  percent  manufacturing), 


best  or  least  shipping  time  from  i  to  j , 

ratio  of  travel  time  by  mode  k  to  the 
best  time. 


best  or  least  cost  of  shipping  from  i  to  j , 

ratio  of  the  cost  of  mode  k  to  the  least 
cost , 


n„  =  number  of  modes  serving  i  and  j  ,  and 
a,B  =  parameters. 
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Baumol,  Studies  on  the  Demand  for  Freight  Trans¬ 
portation  ,  1,  Chapter-!"! 
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The  data  for  this  model  are  disaggregated  by  origin 
and  destination  pairs  and  the  model  contains  a  larger 
range  of  level  of  service  and  market  attributes  than  others. 
However,  its  major  disadvantage  is  the  lack  of  change 
in  mode  preference  due  to  any  improvements  in  performance 
characteristics  of  any  given  mode  short  of  becoming  the 
best  mode  with  that  characteristic .  The  "best”  definition 
may  be  inconsistent  under  some  circumstances,  such  as 
truck  being  the  least  cost  for  very  small  shipments  but 
the  most  expensive  for  very  large  shipments. 

The  Mathematica  study  also  developed  a  model  based 
on  micro  economic  theory  of  the  firm  similar  to  the  general 
inventory  theoretic  cost  function  discussed  in  Chapter  2. 
This  is  a  least  cost  model  based  on  annual  variable  trans¬ 
portation  costs  of  the  industry  in  city  j ,  receiving  com¬ 
modity  q  from  city  i.  Variable  costs  are  defined  as  the 

sum  of  direct  shipper  costs,  total  in-transit  carrying 

22 

costs  and  safety  inventory  costs. 

The  transportation  demand  function  is  as  follows: 


4  8 


where  VC 


annual  variable  transport  related  costs, 
carrying  costs  per  day  (g) , 


shipping  cost  per  unit  of  commodity  q  from 
city  i  to  city  j  by  mode  k, 

time  between  shipments  of  commodity  q  to  j 
the  reciprocal  of  the  number  of  shipments 
per  year , 

quantity  of  commodity  q  shipped  from  city 
i  to  city  j  by  mode  k, 

travel  time  for  shipping  commodity  q  from 
city  i  to  city  j  by  mode  k,  and 

inventory  cost  per  unit  of  commodity  q  per 
year  at  j . 


Kullman  used  a  binary  logit  (logistic)  model  on  aggre¬ 
gate  data  to  predict  the  division  of  traffic  between  truck 
23 

and  rail.  The  model  includes  service,  commodity,  and 
market  attributes  as  follows: 


f  (d^,  P^,  r1/^2,  hi/h2.  al//a2) 


(3.21) 


where  T.  .  .  3 
kij 


d.  . 


quantity  of  commodity  q  moving  by  mode  k 
from  origin  i  to  destination  j 

distance  from  i  to  j , 

value  of  commodity  k  per  ton, 

transport  rate  of  mode  1  or  2, 

transport  time  of  mode  1  or  2,  and 
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Brian  Kullman,  A  Model  of  Rail  and  Truck  Competi¬ 
tion  in  the  Intercity  Freight  Market  (Massachusetts  Insti- 


reliability  of  mode  1  or  2  • 


r 


al ,  2 

This  formulation  is  very  similar  to  that  used  in  the 

discriminant  model  discussed  later  in  Chapters  4,  5  and  6 

and  utilized  in  the  analysis  of  demand  for  transportation 

in  this  dissertation.  Production  economies  of  the  firm 

and  transport  mode  are  included.  The  model  has  been  used 

to  estimate  impacts  of  changes  in  mode  performance  and 

cost,  and  the  impact  of  shifting  relative  market  prices 

between  origin  and  destination. 

Kullman's  model  does  allow  the  evaluation  of  the 

impact  of  changes  in  service  and  rate  attributes  and 

changes  in  reliability,  although  the  aggregate  formulation 

limits  variability  in  the  data  and  reduces  its  empirical 

24 

usefulness  for  estimating  non  aggregated  behavior.  As 


2 >  ^ 

There  is  considerable  literature  reflecting  the  use 
of  logistic  and  probabilistic  (shortened  to  logit  and 
probit)  formulations  of  the  decision-makers  response  to 
choice  of  transport  mode.  Most  of  the  literature  deals 
with  passenger  traffic.  See  P.R.  Stopher,  "A  Probability 
Model  of  Travel  Mode  Choice  to  the  Work  Journey,"  Highway 
Research  Record  283  (1969),  pp .  57-65;  S.  Reichman  and 
P.R.  Stopher ,  "Disaggregate  Stochastic  Models  of  Travel- 
Mode  Choice,"  Highway  Research  Record  309  (1971),  and  Paul 
R.  Rassam,  et  al,  "The  n-Dimens ional  Logit  Model:  Develop¬ 
ment  of  Application,"  Highway  Research  Record  309  (1971), 
pp.  135-147.  The  basic  attributes  of  the  models  are  dis¬ 
cussed  above.  A  more  extensive  discussion  of  the  econo¬ 
metric  characteristics  of  the  model  is  given  in  Chapter  4, 
following. 


stated  in  Chapter  4,  the  linear  logit  model  has  been  criti¬ 
cized  by  Oum  as  inappropriate  to  use  for  transport  demand 
studies  due  to  many  rigid  a  priori  restrictions  on  param¬ 
eters  of  price  responsiveness  of  demand  and  the  underlying 

2  5 

structure  ofpreference  is  irregular  and  inconsistent. 

Discriminant  Models 


Beuthe  developed  a  modal  split  model  based  on  dis¬ 
criminant  analysis  which  accounts  for  the  spatial  distribu¬ 
tion  of  shipments  of  com  from  Illinois  via  barge,  rail 
2  6 

or  truck  modes.  Regional  boundaries  are  defined  by 
the  following  discriminant  model: 

V  -  f  <v  V  v  (3-16) 


where  ^ 

rk 

nk 

<*k 


quantity  of  commodity  q  shipped  by  mode  k 
from  region  j,  where  k  =  1,2,3), 

transport  rate  of  mode  k, 

time  costs  by  mode  k,  and 

quality  differences  between  modes. 
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Tae  Hoon  Oum,  "A  Warning  on  the  Use  of  Linear  Logit 
Models  in  Transport  Mode  Choice  Studies,"  The  Bell  Journal 
of  Economics,  10,  pp.  374-378. 

26 

Michael  V.  Beuthe,  Freight  Transportation  Modal 
Choice:  An  Application  to  Corn  Transportation  (Northwestern 

University,  Evanston,  Ill.,  Dissertation,  1968) ,  Chapter  4, 
pp.  4-78. 
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Beuthe  also  developed  a  spatial  model  by  extension 
of  traditional  location  theory.  Location  of  shipments 
by  each  mode  are  relative  to  distance  and  costs,  and  the 
shipment  predictions  take  the  form  of  boundaries  between 
market  areas  for  the  respective  modes.  The  discriminant 
model  estimated  the  quality  differential  between  modes 
and  provided  substantial  insight  into  the  comparative 
advantage  held  by  each  transport  mode  whereas  the  location 
al  model  was  limited  to  describing  mode  differences  in 
relative  transportation  costs  alone.  The  empirical  study 
was  of  com  shipments  from  Illinois. 

Allen  utilized  discriminant  analysis  to  analyze 

the  division  of  traffic  between  air  and  sea  modes  of  trans 

27 

portation  in  the  North  Atlantic  trade  route.  The  linear 
discriminant  model  was  defined  as  follows: 


(3.17) 


where 


quantity  of  commodity  q  shipped  by  mode  k, 
value  per  pound  of  commodity  q,  and 
rate  of  mode  k. 


Allen  found  that  the  discriminant  approach  worked 


best  on  commodities  and  routes  which  resulted  in  exclusive 
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modal  choice.  Regression  performed  poorly  on  the  entire 

sample.  The  analysis  suffered  because  of  lack  of  data  on 

transit  time  and  shipment  size. 

Antle  and  Haynes  used  linear  discriminant  analysis 

in  a  two  mode  model  for  shipments  in  the  Upper  Ohio  River 
28 

area.  The  model  was  formulated  as  follows: 


where 


s 


kij 


rkij 

ckij 


b  .  . 

aij 


f^kij  *  ^kij  ’  skij  ’  rkij  ’  ckij  ’  baij  ^ 


(3.18) 


?uantity  of  commodity  q  shipped  by  mode  k 
rom  origin  i  to  destination  ] ,  with  k  = 

1  or  2, 

distance  from  i  to  j  by  mode  k, 

size  of  shipment  from  i  to  j  by  mode  k. 

travel  time  from  i  to  j  by  mode  k, 

transport  rate  from  i  to  j  by  mode  k. 

handling  cost  from  i  to  j  by  mode  k,  and 

transport  cost  of  alternative  mode  from 
i  to  j  . 


The  model  was  estimated  for  a  relatively  small  sample 
of  coal,  coke,  petroleum,  and  chemical  shipments  by  barge 
and  rail.  Costs  of  the  alternative  mode  failed  to  appear 
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Lloyd  G.  Antle  and  Richard  W.  Haynes,  An  Applica¬ 
tion  of  Discriminant  Analysis  to  Division  of  Traffic  Be¬ 
tween  Transport  Modes  (U.S.  Army  Engineer  Institute  for 
Water  Resources,  Ft.  Belvoir,  Virginia,  1971),  Report  71-2. 
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in  the  calibrated  model  due  to  lack  of  significance.  Th 
two  group  procedure  developed  by  Tintner  was  utilized  as 
follows : ^ 

Z  =  k,  d-,  +  kod0  4-  ....  +  k  d  (3.1 

11  L  L  p  n 

where  Z  =  test  value, 

kp  =  linear  weights,  and 

d  =  differences  in  means  in  each  group  of  n 
variable. 

The  weights  (k^'s)  are  assigned  in  a  way  to  maximize  the 

2 

value  of  Z  relative  to  its  variance.  The  resulting 


p(x) 


u  X  (weighted  vector  of 

characteristics) 

Figure  3.1  Probability  of  occurrence  in  groups 
1  or  2  depending  on  weighted  charac¬ 
teristics  &  Discriminant  Function  (Z) 

^^Gerherd  Tintner,  Econometrics . (John  Wiley  &  Sons 
1964) ,  Chapter  6 . 


computations  are  similar  to  linear  regression  and  classifi¬ 
cation  is  based  on  a  critical  value  of  Z  which  minimizes 
classification  errors. 

The  resulting  model  performed  relatively  well  pre¬ 
dicting  the  demand  for  coal  shipments,  but  the  sample 
was  too  small  to  permit  good  prediction  for  other  commodity 
groups.  The  discriminant  function  allowed  an  estimate 
of  the  demand  function  by  means  of  parametric  shifting 
of  barge  rates  to  derive  the  quantity  shipped  by  barge 
to  zero.  The  model  allows  an  assessment  of  the  impacts 
of  rate  change  and  other  attributes  on  the  choice  of  mode 
and  therefore  the  demand  by  mode. 

Southwestern  Division,  Corps  of  Engineers,  applied 

discriminant  analysis  to  selected  commodity  shipments  on 

30 

the  Arkansas  River  area.  The  Southwestern  Division 
analysis  was  based  on  194  shipments  data  by  rail,  truck 
and  barge  modes  for  six  commodity  groups  which  were  moving 
by  barge  transportation  for  the  year  1971.  The  analysis 
found,  essentially,  a  highly  price  inelastic  demand  func¬ 
tion  for  water  transportation  based  on  the  data.  Since 
the  data  were  collected  in  the  first  year  that  navigation 

30 

Southwestern  Division,  Discriminant  Analysis 
Applied  to  Commodity  Shipments  in  the  Arkansas  River  Area , 
Report  74-R2  (U.S.  Army  Engineer  Institute  for  Water 
Resources,  Ft.  Belvoir,  Virginia,  1974). 


was  fully  operational  on  the  Arkansas  River,  one  of  the 
possible  explanations  for  the  inelastic  properties  of 
the  demand  function  is  that  the  shippers  were  attempting 
to  stimulate  downward  reductions  in  rail  rates  and  that 
was  a  prominent  part  of  the  motivation  for  using  water. 

This  hypothesis  was  confirmed  later  in  a  study  of  rail 
rate  adjustments  which  showed  that,  at  least  in  iron  and 
steel  products,  substantial  drops  in  rail  rates  were  begun 
in  1971.  The  model  adopted  was  exactly  the  same  as  was 
used  in  the  study  by  Antle  and  Haynes  in  1971. 

Sasaki  developed  a  comprehensive  commodity  demand, 

supply  and  transport  demand  model  to  project  the  demand 

31 

for  freight  transport  by  mode.  The  comprehensive  model 
was  composed  of  three  submodels.  Model  1  projects  output 
of  firms  which  utilize  coal  as  a  factor  of  production. 

Model  2  allocates  the  available  supply  of  coal  in  a  way 
that  minimizes  transportation  costs  to  satisfy  demands. 
Model  3  is  a  mode  choice  model  which  introduces  nonmone- 
tized  mode  choice  variables . 

The  discriminant  mode  choice  model  was  developed 
for  rail,  unit  train  and  barge  modes  as  follows,  with 
the  variables  defined  as  in  3.18: 

31 

Brion  R.  Sasaki,  A  Regional  Model  of  the  Future 
Demand  for  TransportationT  The  Case  of  Barge  Transporta¬ 
tion  (University  of  Cincinnati,  Cincinnati,  Ohio,  Disserta- 
tion,  1974) . 
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(3.20) 


T,  .  .  ^  =  f  ( d,  .  .  ,  h,  .  .  ,  s,  .  .  ,  c,  .  . ) 
kij  kij  kij  ’  kij  ’  kij 

Sasaki's  model  was  estimated  from  97  observations 
in  three  Bureau  of  Economic  Analysis  (Department  of  Com¬ 
merce)  economic  regions  in  the  Ohio  Valley.  His  work 
confronts  the  primary  problems  in  projecting  freight  de¬ 
mands;  changing  technology  of  product  users  (now  compli¬ 
cated  by  environmental  regulations) ,  supply  constraints 
which  affect  the  quantity  available  (also  affected  by 
environmental  regulations) ,  and  shifting  service  and  cost 
characteristics  of  transport  modes,  but  does  not  deal 
with  long-run  locational  changes  which  would  affect  trans¬ 
port  demand. 

The  comprehensive  model  contains  some  redundancy 
in  the  allocation  of  aggregate  shipments  from  origin  to 
destination  in  model  2  and  the  mode  choice  calculation 
in  model  3.  The  allocation  and  choice  of  mode  could  have 
been  made  in  one  step  if  all  variables  had  been  converted 
to  total  cost  in  advance  of  the  cost  minimization  calcula¬ 
tion.  This  suggests  a  possible  use  of  modal  split  in 
an  earlier  step  to  determine  the  implicit  prices  that 
would  be  assigned  to  nonmonetary  variables  in  selection 
of  mode.  Shipment  size,  however,  would  resist  monetization 
but  perhaps  is  reflected  in  the  scale  economy  of  rates. 
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The  model  is  capable  of  evaluating  the  impact  of 
changes  in  service  and  rate  attributes  and  is  uniquely 
able  to  accept  constraints  in  both  the  demand  and  supply 
functions  for  coal. 

Comparative  Studies 

Hartwig  and  Linton  developed  logit,  probit  and  dis¬ 
criminant  functions  to  model  individual  shipper  choice 

of  mode  between  full  load  truck  and  full  load  rail  for 

32 

intercity  freight  movements.  The  model  was  estimated 
from  1213  freight  way  bills  for  truck  and  rail  shipments 
of  consumer  durables  as  follows : 

f  (Pq.  rUj  -  r2ij  ’  sl  -  s2'  alij  -  a2ij> 

(3.22) 

quantity  transported  by  rail  (1)  or 
truck  (2)  from  origin  i  to  destination  j, 


rl .  2ij 

S1 , 2 
al,2ij 


price  per  ton  of  commodity  q, 

actual  transport  rate  by  mode  l-,  2  from 
i  to  j, 

actual  shipment  size  by  mode  1,  2,  and 
reliability  of  mode  1,  2  from  i  to  j . 


J.  Hartwig  and  W.  Linton,  Disaggregate  Mode  Choice 
Models  of  Intercity  Freight  Movement  (Nortnwestem  Univer- 
sity,  Evanston,  Illinois,  Thesis,  1974) . 
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They  found  that  the  probit  and  logit  formulation 
performed  somewhat  better  than  discriminant  functions 


for  the  commodities  selected.  The  model  allows  for 
evaluation  of  impacts  from  shifts  in  transport  costs  and 
reliability  measures  but  does  not  include  commodity 
attributes  or  length  of  time  in  transit  measures. 

Summary 

Analysis  of  the  literature  in  freight  demand  modeling 
shows  the  steady  progress  from  dependence  on  aggregate 
data  models  which  tended  to  reduce  applicability  to  par¬ 
ticular  demand  and  suppresses  variability,  to  models  using 
a  much  larger  set  of  variables  which  influence  the  amount 
of  freight  to  be  shipped,  its  origin  and  destination,  and 
the  choice  of  transport  mode.  Models,  furthermore,  which 
concentrate  on  traffic  flows  and  transport  rates  miss 
important  explanatory  variables.  Behavioral  models 
in  which  individual  shipper  decisions  are  optimizing 
choices  are  more  likely  to  produce  informed  analysis  of 
the  causal  variables  relating  to  quantities  shipped. 

Mode  choice  models  (whether  estimated  by  some  form  of 
regression  or  by  discriminant  analysis)  offer  the  possi¬ 
bility  of  estimating  both  the  quantity  of  traffic  to  be 


shipped  by  a  mode  and  the  demand  function  for  that  mode 
(that  relationship  between  transport  rate  and  quantity 
shipped) .  A  description  of  these  models  is  in  the  next 
chapter. 
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CHAPTER  IV 


ECONOMETRIC  CHARACTERISTICS  OF 
REGRESSION  AND  DISCRIMINANT  METHODS 


Introduction 


The  following  chapter  discusses,  in  turn,  the  econo¬ 
metric  characteristics  of  several  varieties  of  regression 
and  discriminant  analysis.  Linear  regression  is  used  to 
fit  the  parameters  of  gravity,  aggregate  demand  and  ab¬ 
stract  mode  models.  Nonlinear  regression  models  are 
utilized  in  probit  and  logit  form  of  modal  split  models. 

It  is  shown  that  the  linear  discriminant  analysis  of  a 
two  group  case  is  a  special  case  of  linear  regression. 

Linear  Regression 

Following  Johnston,  assume  a  linear  relationship 
between  variable  Y  and  k-1  explanatory  variables  X2 ,  X^ 
through  X^  and  a  disturbance  term  u.^  With  a  sample  of 
n  observations  on  y  and  x's  the  following  can  be  written: 


Yi  " 


bi  +  b2X2i 


Vki  +  ui 


(4.1) 


^•J.  Johnston,  Econometric  Methods  (McGraw  Hill,  1960) 
See  Chapters  1  and  2  for  the  two  variable  linear  model 
and  Chapter  4  for  the  general  linear  model. 
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The  coefficients  and  disturbance  are  unknown  and 
the  problem  is  to  obtain  estimates  of  the  unknown.  In 
matrix  notation  the  results  are  the  following: 


Y=X8+u  (4.2) 
The  simplest  set  of  assumptions  are: 

E(u)  =  0  (4.3a) 
E(uu')  =  o2ln  (4.3b) 
X  is  a  set  of  fixed  numbers  (4.3c) 
X  has  rank  k  <  n  (4.3d) 


It  is  assumed  that  the  expected  value  of  the  u. ’s  is 
zero,  that  the  u^ 1 s  have  constant  variance  and  that  off 
diagonal  elements  are  pairwise  uncorrelated,  that  the  sole 
source  of  variation  on  the  Y  vector  is  variations  in  the 
U  vector  and  that  the  properties  of  the  estimators  and 
tests  are  conditional  upon  the  X’ s  and  that  the  number  of 
observations  exceeds  the  number  of  parameters  to  be  esti¬ 
mated  and  no  exact  linear  relationship  exists  between  the 
X  variables. 

Least  square  estimates  of  the  form: 

Y  =XB+  e  (4.4) 

can  be  derived  by  minimizing  the  square  of  the  e^  term, 

t 
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(4.4a) 


e*  =  ('<  =  (Y-XB)'  (Y-XB) 


=  Y'Y  -2  B'X'Y-  B'X  XB 

(4.4b) 

then  finding  the  value  of  B  which  minimizes 

uals  by  taking  the  derivative: 

squared  resid- 

(e'e)s-2X'Y+2X'XB 

(4.5) 

then  X#XB  =  X#Y  . 

(4.6) 

and  B  =(X/X)*1(X/Y)  . 

(4.7) 

The  least  square  estimates  are  best  linear  unbiased 

In  the  special  case  where  the  joint  distribution  of  X  and 

Y  is  multivariant  normal,  the  least  square  estimate  has 

the  properties  of  consistency,  efficiency,  minimum  vari- 

3 

ance,  unbiased,  and  sufficiency. 

Estimation  by  maximum  likelihood  methods  can  accom¬ 
modate  nonlinear  relationships  in  the  parameters  but 
require  assumptions  of  normality  of  the  u^ ' s  : 

u-_  -  N  (0,  a2  In)  (4.8) 


Form  the  likehood  function: 


L 


_ 1 

(2tto2 


)U/2 


(4.9) 


2Ibid. ,  p.  111. 

Ibid. ,  p.  133,  and  Franklin  A.  Graybill,  An  Intro¬ 
duction  to  Linear  Statistical  Models  (McGraw-HilTi  1961 ) ; 
Vol .  1,  Chapters  5- 10 . 
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(4.10) 


r 


exp 


_  (Y-XB)*  (Y-XB) 

L  2a2  J 


||  =  (2d2)'1  2X*  (Y-XB)  =  0 


.2,-1 


(4.11) 


|^2  =  §  (o2)'1  +  \  (a2)'2  (Y-XB)1  (Y-XB)  =  0, 

(4.12) 

simplify  and  solve  for  the  maximum  likelihood  estimators 

A  A 

B  and  o  -. 

B  =  (X'X)1  (X*Y),  (4.13) 

~2  =  n'^Y-XBj'tY-XB).  (4.1.4) 

Any  estimate  is  equal  to  plus  a  linear  function 
of  u,  which  has  a  multivariate  normal  distribution,  there¬ 
fore  B^  has  a  normal  distribution.  This  allows  the 
derivation  of  significance  tests  and  confidence  intervals 
for  B^.  The  critical  assumption  for  regression  analysis 
is  that  the  error  term  is  distributed  normally.  Procedures 
have  been  developed  to  cope  with  problems  related  to 
heteroskedasticity  and  the  variety  of  cases  where  the  inde¬ 
pendent  variables  are  not  pairwise  uncorrelated. 


Logistic  Regression  (Logit) 

Logit  analysis  is  based  on  the  premise  (from  studies 
of  toxicology)  that  members  of  a  population  subjected  to 
a  stimulus  which  can  range  over  an  infinite  scale  will 
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respond  (in  a  binary  choice  situation)  according  to  a 
simple  sigmoid  curve. ^  In  the  context  of  choice  or  trans¬ 
port  mode,  the  division  of  traffic  between  two  modes  is 
between  0  and  1  and  their  sum  is  unity.  Mode  choice  is 
a  monotonic  function  of  independent  variables  and  the 
transportation  variables  are  expressed  in  units  such  that 
if  there  is  an  increasing  or  decreasing  utility  of  trans¬ 
portation  by  a  given  mode,  then  the  share  of  that  mode 
decreases  or  increases  when  any  of  its  transportation 
variables  increases  or  decreases  relative  to  the  other 
mode.  Figure  4.1  shows  the  assumed  response  function. 


PERCENT  MARKET 

SHARE  SHIPPED  BY  MODE  i 

100  1 - 


STIMULUS 

Figure  4.1  Cumulative  Response  to 
a  Varying  Stimulus. 


^Reichman  and  Stopher,  "Disaggregate  Stochastic 
Models  of  Travel-Mode  Choice,"  p.  95. 
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This  nonlinear  response  can  be  written: 


Pi 


1-e  §(x) 


(4.14) 


where  P . 

i 


the  conditional  probability  of  choosing 
mode  i ,  and 


g(x)  =  a  utility  function. 


In  choosing  between  two  transport  modes,  it  is 
assumed  in  the  logit  model  that  the  reaction  to  stimulus 
is  determined  by  the  utilities  of  each  mode  and  the  charac¬ 
teristics  of  the  user.  In  the  two  mode  case  the  logit 
probability  function  and  the  resultant  logit  model  can  be 
written  as  follows:^ 


P  ( i) 


expressed  as, 


/  N  \ 

exp  ^-a0-a.x.-2anxnj 


1  + 


,  which  can  be 

(4.15) 


108  (l^FTT))  '  ao  +  aixi  +  an  V 
'  '  n -z 


(4.16) 


where  P(i) 


x . 

l 


n 


n 


=  the  conditional  probability  of  choosing 
mode  i, 

=  the  price  of  mode  i  relative  to  the  other 
modes , 

=  other  exogenous  factors  affecting  mode 
choice  n=2,3...N,  and 

=  the  parameter  of  the  logistic  function 
n  -  0,  1  ...  N. 


^Oum,  "A  Warning  on  the  Use  of  Linear  Logit  Models,” 
p.  375. 


The  parameters  are  estimated  by  linear  regression. 

The  dependent  variable  is  the  logarithm  of  the  market  share 
ratio  of  two  modes.  Two  forms  of  the  utility  function 
have  been  postulated.  Ranson  et  al,  and  Kullman,  specified 
the  price  variable  as  a  ratio  of  the  price  of  the  ith  mode 
to  the  base  mode.^  McFadden  along  with  Richards  and  Ben- 
Akiva  expressed  the  price  variable  as  the  difference  be¬ 
tween  the  price  of  the  ith  mode  and  the  base  mode.^  If 
the  response  function  shown  in  Figure  4.1  represents  the 
cumulative  normal  distribution,  the  resulting  probability 
(probit)  model  is  estimated  by  maximum  likelihood  methods. 

Discriminant  Analysis 


The  discussion  begins  with  the  two  group  case  pre- 

Q 

sented  by  Morrison,  then  moves  to  the  k  group  case  pre- 

Q 

sented  by  Tatsuoka.  The  function  of  discriminant  analysis 
is  to  form  a  linear  combination  of  predictor  variables 


Rassam,  et  al,  "The  N-Dimensional  Logit  Model: 
Development  and  Application";  Kullman,  A  Model  of  Rail/ 
Truck  Competition  in  the  Intercity  Truck  Market ;  and 
others . 

^D.A.  McFadden,  "The  Measurement  of  Urban  Travel 
Demand,"  Journal  of  Public  Economics,  Vol.  3  (1974),  pp . 
303-328;  and  M.G.  Richards  and  Moshe  E.  Ben-Akiva,  A  Dis¬ 
aggregated  Travel  Demand  Model  (Lexington  Books,  1975). 

Q 

Donald  G.  Morrison,  "On  the  Interpretation  of  Dis¬ 
criminant  Analysis,"  Journal  of  Marketing  Research,  b 
(May  1966);  pp .  156-163. 

q 

Maurice  M.  Tatsouoka,  Multivariate  Analysis:  Tech¬ 
niques  for  Educational  and  Psychological  Research  (John 


Figure  4.2  Regression  and  Discriminant  Functions, 

(a)  Regression  Where  Group  Membership 
is  Indicated  By  a  Dummy  (1,0)  Variable 

(b)  Discrimanant  Function  for  a  Two 
Group  Case. 


Figure  4.3  Discriminant  Functions  for 
the  Three  Group  Case 

In  Figure  4.2  and  4.3  the  cocrdinants  represent 
variables  which  distinguish  between  groups.  In  the  trans 
port  case  each  group  represents  a  transport  mode  and  one 
variable  could  be  costs  and  the  other  service.  With  this 
model  of  mode  choice,  discriminant  analysis  derives  the 
boundaries  which  minimize  the  group  members  which  would 
be  classified  in  error. 
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For  notation  let, 


x 


ji 


the  ith  group  members  value  of  the  j th 
independent  variable, 


bj 


the  discriminant  coefficient  for  the  j th 
variable , 


Z^  =  the  ith  group  members  individual  discrimi¬ 

nant  score,  and 


Z 


crit 


critical  values  of  the  discriminant  func¬ 
tion. 


Then  let  each  group  member's  discriminant  score  Z ,  be 
a  linear  function  of  the  independent  predictor  variables, 


=  b. 


+  Vu 


+  b2*21 


b  x  .  . 
n  ni 


(4.17) 


If  Zi  >  Zcrit’  classify  1  Co  group  1  and  if  <  Zcrit, 
classify  i  to  group  2.  The  locus  of  points  where, 


+  vu 


,  +  b  x  . 
n  ni 


=  Z 


crit  ’ 


(4.18) 


would  be  the  classification  boundary.  In  the  two  group 
case  the  classification  boundary  is  a  two  dimensional 
plane  in  three  dimensional  space.  In  general  the  classifi¬ 
cation  boundary  is  an  n-1  dimensional  hyperplane  in  n 
space . 

The  following  discussion  develops  the  general  dis¬ 
criminant  analysis  case,  continuing  with  the  explanation 
by  Tatsuoka.  The  first  step  is  to  decide  on  a  criterion 
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for  measuring  group  mean  differences.  Th:  F  ratio  could 


serve  as  an  appropriate  criterion.  For  K  groups  with  a 
total  of  N  individuals,  the  F  ratio  is  given  by: 


F 


SSb/(N-l) 
SS  /(N-K) 

W 


SSb  N-K 

SS  K-l 
w 


(4.19) 


where  =  sum-of -squares  between  group  means, 

and  SSw  -  suic-of- squares  within  group  means. 

N-K 

The  second  term  (  ^k-1^  a  constant  for  any 

given  problem,  therefore,  the  essential  term  is  the  rela¬ 
tion  SS,  /SS  . 

D  W 

With  p  predictor  variables,  ,  X2  •  ••  X  ,  form  the 
linear  combination: 


Y 


V1X1  +  V9X2 


v  X 
P  P 


(4.20) 


where  the  v.'s  are  linear  weights,  and  for  which  the  within 
group  means  and  between  group  means  sum-of- squares  are  ex- 
pressable  as  quadratic  terms. 

Denote  the  sum  of  squares  of  Y  within  the  kth  group 
by  SS^  (Y)  and  let  v’  =  (v2,v2  ...  v  ) ,  then, 


SSw(Y)  =  SS1(Y)  +  SS2(Y)  +  ...  SSk(Y) 


where  SSk<Y) 


S,  S, 

k , yy  k , yx 

S,  S, 

k , xy  k , xx 


*k 


(4.21) 


(4.22) 
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Then  SS(Y)  = 
w 

J2W+.  .  •  ¥'Sk¥ 

= 

w/(S1+S2  +  -  •  ■  Sk)v 

(4.23) 

or  SS  ( Y) 
w 

v'Nv 

(4.24) 

k 

since  VS.  = 

W,  by  definition  in  (4.21). 

(4.25) 

II 

k=1 

To  obtain 

the  formulation  for  between-group 

means 

sum-of -squares ,  define  the  diagonal  elements  of  B  as  the 
usual  between-group  means  sum-of -  squares  for  the  variables 
taken  one  at  a  time: 


=  I  nk  (Xik-X^2  for  all  i, 
k  1 


where  nk 


=  size  of  kth  group, 

=  kth  group  mean  of  X^, 
=  grand  means  of  X^ . 


(4.2G) 


The  off  diagonal  elements  of  B  are  then  between  group  mean 
sums  of  products  for  pairs  of  variables  over  the  K  groups: 

bij  ■  nk  <-\k  -  V  (Xjk  -  V  fori#b 

(4.27) 

Changing  the  notation  to  a  more  general  form  let 


X  =  group  means  and, 

X  =  grand  means , 
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chert,  B 


(4.28) 


=  (X-X)'(X-X) 


Pre  and  post  multiply  both  sides  of  (4.28)  by  v*  and  V  re¬ 
spectively  , 

v'Bv  =  v'(X-X)'  (X-|)v  (4.29) 

=  (Xv  -Iv)'  (Xv-Xv)  (4.30) 

Since  Xv  and  Xv  are  vectors  which  elements  consist  of 

the  group  means  of  Y  and  the  grand  means  of  Y  respectively , 

the  product  is: 

(Xv-Iv)'(Xv-lv)  =  l  n,  (Y,  -  Y)2  (4.31) 

k=l 

where  Y,  =  group  means  of  Y, 

and  Y  =  grand  mean  of  Y, 

which  is  the  between-group  means  sum-of-squares  of  the 
transformed  variable  Y.  Thus, 

SSb(Y)  =  v'Bv  (4.32) 

and  rewriting  (4.24)  and  (4.32)  gives; 


SSb(Y) 

ssTTyT 

w 


v'Bv 

v'Wv 


(4.33) 


The  ratio  A  (defined  in  (4.39))  provides  a  criterion  for 
measuring  group  differences  in  the  dimensions  specified 
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by  vector  v.  This  is  the  discriminant  criterion  proposed 
by  Fisher  in  1936.^ 

The  next  step  is  to  maximize  the  discriminant  cri¬ 
terion  (subject  to  v'v  =  1.0)  by  determining  a  set  of 
weights  (v^,  ...  vn)  which, 

3X  _  2f(Bv)(w,Wv)-  \(v*Bv)(WvVl  =  0  (4.34) 

3V  Iw'WwT^ 

Simplifying  the  equation  reduces  to: 

'  ***f^f4.  ,  o.  (4.35) 

which  implies 


Bw  -  xWw  =  0, 
since  (v'Nv)  =  0  . 

(B  nwlv  =  0. 


(4.36) 


(4.37) 


Assume  W  is  nonsingular  and  therefore  possesses  an  inverse, 
premultiply  both  sides  by  W 

( W"1  B- A  ■ )  v  =0.  (4.38) 


The  solution  yields  eigenvalues  A^  and  associated  eigen¬ 
vectors  w  which  maximize  the  discriminant  criterion, 
m 


^R.A.  Fisher,  "The  Use  of  Multiple  Measurements  in 
Taxonomic  Problems,"  Annuals  of  Eugenics,  Vol.  7  (1936), 
pp.  179-188. 
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The  eigenvector,*^  =  (v-q  ,  v-^  •••  )  correspond¬ 

ing  to  the  largest  eigenvalue,  A^,  provides  a  set  of 
weights  such  that  the  transformed  variable, 


V11X1  +  V12X2 


+  v-t  X 
lp  p 


(4.39) 


has  the  largest  discriminant  criterion  A-^  attainable  by 
any  linear  combination  of  p  predictor  variables. 

In  the  two  group  case,  computation  of  the  single 
discriminant  function  is  considerably  simplified.  In 
particular  B  assumes  a  particularly  simple  form.  The 
expression  reduces  to: 


bij  *  nl  +n2 

(4.40) 

and  the  grand  mean  X  is  equal  to: 


nl  Xil  +  n2  xi2 
n^  +  n2 


(4.41) 


and  substituting: 

nl  n2  - 

b,  .  =  z~~  (X.,  -  X,0)  (X.,  -  X.0l 
ij  n^+n2  *  tl  i2;  v  jl  j2J 

Define  a  p  dimension  row  vector, 


(4.42) 


^  [*11  "  *12  ’  *21  '  *22  ••••  *pl 


Xp2]  (4.43) 


of  the  differences  between  the  group  means  of  the  p 
variates,  then 
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n1n? 

B  =  — L-L.  dd' 

nl  +  n2 


(4.44) 


and  define  c  = 


nln2 
nl  +  P2 


,  substitute  into  4.38, 


[cW_1(dd')-ul  J  w  = 


(4.45) 


where  u  =  since  equation  4.38  can  be  rewritten: 


Bv  =  A  Wv, 


(4.46) 


then  add  ABv  to  both  sides; 


(l+A)Bv  =  A  (W+  B)  v  ,  or 
(T'’B  -  ITT  «)«■  0. 


(4.47) 

(4.48) 


where  T  is  the  total  sum  of  squares  and  cross  products 
matrix . 

Tatsuoka  shows  that  if  V0  is  the  eigenvector  of 
W-1B  associated  with  the  eigenvalue  A  ,  it  is  also  the 
eigenvector  ofT”^B,  and  its  associated  eigenvalue  is 
AQ/1  +  AQ. 1L  Thus  the  eigenvalues  stand  in  relation; 

A 


u 


1  +  A 


(4.49) 


11 


Tatsuoka,  Multivariate  Analysis,  pp .  172-174. 
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uv 


(4.50) 


=  c  T'1  (dd')  v  ( 

=  c  T"1  (d)  (d'v)  .  (4.51) 

Since  d*v  is  a  scalar,  but  unknown  quantity  (d*  it 
lxp  and  y  =  p  xl) ,  collect  all  scalar  quantities  into  a 
single  multiplier  and  write: 


v  =  [c(d/v)/u]T'1d,  and  (4.52) 

v  =  mT^d  (4.53) 


where  y  is  an  unknown  scalar,  since  it  is  pxl  and  m  = 

[  (d'v)]/u  .  This  gives  a  solution  for  the  eigenvector 
y  o£T*^B  in  the  two  group  case,  because  v  is,  in  any  case, 
proportional  up  to  an  arbitrary  proportionality  constant. 

In  the  two  group  case,  the  single  discriminant  function 
can  be  solved  without  solving  the  eigenvalue  problem  by 
postmult  ip lying  the  inverse  of  the  sums  of  squares  and 
cross  products  (SSCP)  matrix  by  the  column  of  mean  differ¬ 
ences  on  the  p  predictors  to  get  the  vector  of  discriminant 
coefficients.  The  single  discriminant  function-  can  be 
obtained  without  solving  the  eigenvalue  by  normalizing 
y .  Rewrite (4 . 53) , 


v 


S'1  md 

PP 


(4.54) 


where  S  equals  the  sum  of  squares  and  cross  products 
of  other  predictor  variables.  If  the  criterion  variable 
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—ha 


4 


m 


in  the  cwo  group  case  is  a  dichotomous  variable  (1,0), 
we  can  infer  that  the  sum  of  the  products  is  given  by: 

2>y  =  nln2  (X,  -  X2)  (4.55) 

nl+n2 

Thus  the  elements  of  mdf  given  in  4.54  are  proportional 
to  the  predicted  criterion  sum  of  products  where  the  cri¬ 
terion  dichotomous  variable  indicates  group  membership. 
Therefore  in  the  two  group  case  the  discriminant  weights 
are  proportional  to  the  weights  of  a  multiple  regression 
equation  of  a  dichotomous  group  membership  variable  on 
the  p  predictors.  Several  writers  have  stated  that  dis¬ 
criminant  analysis  in  general  is  a  special  case  of  multiple 
regression,  but  the  statement  is  true  only  for  the  two 
group  case.  In  cases  of  greater  than  two  groups,  dis¬ 
criminant  analysis  reduces  to  canonical  correlation  analy- 
13 

sis.  This  is  possible  if  the  criterion  variables  are 
dummy  variables  with  one  fewer  dummy  variables  than  the 
number  of  groups.  Then  the  predictor  and  criterion  sets 
are  treated  by  the  means  of  canonical  analysis. 

1 2 

K.E.  Garret,  "The  Discriminant  Function  and  Its 
Use  in  Psychology,"  Psychometrics ,  Vol.  8  (1943),  pp .  65- 
79,  and  G.S.  Wherry,  "Multiple  Bi-serial  and  Mutiple  Point 
Bi-Serial  Correlation,"  Psychometrics,  Vol.  12  (1947), 
pp.  189-195. 

13 Ibid. ,  p.  173. 
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Canonical  Analysis 


Canonical  analysis  is  a  technique  to  determine  a 
linear  combination  of  p  predictors  and  a  linear  combination 
of  q  criterion  variables  such  that  the  correlation  between 
these  linear  combinations  in  the  total  sample  is  as  large 
as  possible.  The  problem  is  to  determine  one  set  of 
weights,  u'  =  ( u-j_  ^  ...  u^)  for  the  predictor  variables 
and  another  set  of  weights  V*  =  (v-^,  v ^  •••  v  )  for  the 
criterion  variables  in  such  a  way  the  correlation  r  be- 
tween 


Z 

-  Uq  X-^  +  + 

• • •  +  X^ 

P  P 

(4, 

.56) 

and 

W 

=  V1  V  v2Y2  + 

. . .  X  Yn 

q  q 

(4. 

,57) 

is  the  largest  attainable  for  the  sample.  Partition  the 
sums  of  squares  and  cross  products  matrix  (S) 


s  « 

s 

pp 

'pc 

s 

c 

L  cp 

'ccj 

(4.58) 


then  compute  the  quadruple  matrix  product: 

A=  S'1  S  S"1  S  (4.59) 

#pp  pc  cc  cp 

The  eigenvalue's  u^  and  the  eigenvectors  of  the 

2 

matrix  A  are  computed.  The  largest  eigenvalue  u^  is  the 
square  of  the  maximum  r2w  and  is  the  first  canonical 
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correlation  between  predictor  and  criterion  sets.  The 


elements  of  the  corresponding  eigenvector  Vj  are  the 
weights  to  be  used  in  combining  the  predictor  variables 
to  obtain  the  optimum  linear  combination. 

The  resulting  linear  correlation  Z^,  Z2  ...  Z^ 
are  identical  (within  proportionality)  to  the  discriminant 
function  obtained  in  maximizing  the  discriminant  criterion 
Tatsuoka  proved  that  the  discriminant  criterion  value 

X  is  related  to  the  corresponding  squared  canonical  corre- 

2  14 

late  value  u^  by  the  eigenvector: 

\L  =  ui2/(l-ui)2  (4.59) 

Summary  and  Preview 

Having  reviewed  the  theory  of  transport  demand,  the 
various  approaches  to  demand  modeling  and  the  econometric 
assumptions  of  regression  analysis,  discriminant  analysis 
and  canonical  correlation,  the  choice  was  made  of  the  tech 
nique  to  analyze  transport  demand  for  this  dissertation. 
The  discriminant  analysis  technique  was  chosen  because 
discriminant  analysis  can  provide  a  modal  choice  which 
can  simulate  the  demand  schedule  for  any  given  mode  in 
the  model.  The  schedule  is  estimated  without  any  a  priori 
restrictions,  common  to  linear  regression.  Discriminant 

14 

Maurice  M.  Tatsuoka,  The  Relationship  Between 
Canonical  Correlation  and  Discriminant  Analysis,  Cambridge 
Massachusetts,  Educational  Research  Corporation  (1953). 


analysis  avoids  the  difficulties  embodied  in  linear  logic 
models  due  to  the  rigid  a  priori  restrictions  on  price 
responsiveness  of  demand  and  due  to  irregularity  and  in¬ 
consistency  of  the  structure  of  technology  and  preferences 
underlying  the  linear  logit  model.  On  the  other  hand,  dis¬ 
criminant  analysis  requires  the  assumption  that  each  group 
is  equal  size,  and  is  distributed  normally  with  equal 
variance-covariance  matrices.  Linear  regression  requires 
the  assumption  of  a  normal  distribution  of  the  error  term, 
homoskedasticity  and  independence  of  the  predictor  vari¬ 
ables.  Certain  options  are  available  to  cope  with  situa¬ 
tions  in  which  they  are  not  met  by  the  data.  The  data 
utilized  in  this  analysis  proved  to  be  skewed  but  trans¬ 
formed  comfortably  into  lognormal  form.  Adjustment  to 
correct  the  classification  bias  inherent  in  unequal  size 
groups  is  available  and  was  utilized. 

Applications  of  discriminant  analysis  to  be  applied 
economic  analysis  has  increased  sharply  in  the  past  decade. 
The  ability  to  classify  credit  risks,  the  identification 
of  completed  income  tax  forms  to  be  audited  and  application 
to  mode  choice  models  are  some  examples.^ 

1 5 

See  Edward  I.  Altman,  "Financial  Ratios,  Discrimi¬ 
nant  Analysis  and  the  Predictor  of  Corporate  Bankruptcy," 
The  Journal  of  Finance  (1968)  23;  pp.  589-609;  and  Morrison, 
"On  the  Interpretation  of  Discriminant  Analysis,"  p.  156. 
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The  next  chapter  shows  the  results  of  analysis  of 
several  commodities  moving  by  water  and  competing  transport 
modes  in  the  Ohio  River  and  Arkansas  River  areas . 
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CHAPTER  V 


APPLICATION  OF  DISCRIMINANT  ANALYSIS  MODEL 
OF  TRANSPORT  DEMAND 

This  chapter  discusses  application  of  discriminant 
analysis  to  estimate  demand  functions  for  water  transporta¬ 
tion  in  the  Ohio  and  Arkansas  River  area.  The  model  is 
calibrated  on  a  specially  altered  discriminant  model  which 
first  fits  the  discriminant  functions.  Then  it  proceeds 
to  simulate  the  transportation  demand  function  by  incre¬ 
menting  (adding)  transport  costs  (or  any  other  predictor 
variable  such  as  time  in  transit)  for  any  mode  specified. 
The  original  program  was  the  BMD07M  version  from  the  Bio¬ 
medical  (BMD)  package  developed  by  the  University  of 
California,  Berkley.  The  Statistical  Programs  for  Social 
Sciences  (SPSS)  ,  which  gives  similar  results  except  for 
the  simulation  feature  was  also  utilized.  Both  programs 
are  capable  of  computing  transformations  of  variables. 

Each  develops  a  statistical  summary  of  the  data,  fits 
the  discriminant  functions,  classifies  each  observation 
and  displays  a  summary  of  the  classification  results. 

The  modified  BMD07M  program  performs  the  incrementing 
step  by  adding  prespecified  increases  in  rates  (or  other 
variables),  and  reclassifies  observations  based  or  the 
discriminant  functions  fit  on  the  original  data.  This 


33 


allows  the  simulation  of  the  demand  function  for  any  trans¬ 
port  mode  represented  in  the  data.  The  following  sections 
discuss  the  results  obtained  in  the  analysis  of  data  on 
commodities  which  are  shipped  by  water  and  competing  trans¬ 
port  modes . 

Demand  Model 
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discussion  begins  with  aggregated  coal  shipments  disaggre¬ 
gate  by  type  of  user  and  destination,  then  moves  to  other 
commodities.  A  detailed  summary  of  the  data  used  in  this 
analysis  is  shown  in  Appendix  A. 

Aggregate  Level.  Table  5.1  shows  the  data  from  a 
total  of  395  coal  shipments.  Five  modes  transported  annual 
tonnage  in  excess  of  122  million  tons.  Barge  and  rail 
barge  modes  haul  for  61  percent  of  the  traffic  in  the 
sample  while  regular  and  unit  trains  haul  for  38  percent, 
leaving  about  one  percent  of  the  traffic  to  be  hauled  by 
truck. 

TABLE  5.1 

Summary  of  Coal  Traffic  by  Transport  Mode 
Ohio  River  and  Arkansas  River  Areas 


Rail 

Barge 

Truck 

Rail 

Barge 

Unit 

Train 

No.  of  Shipments 

257 

68 

26 

23 

18 

Av.  Annual  Shipments 
(1000  Tons) 

81 

1,013 

27 

263 

1,418 

Avg  Shipments  (Tons) 

361 

3,831 

50 

1,590 

2,670 

Length  of  Haul  (Miles) 

157 

118 

51 

308 

67 

Time  of  Transit  (Hours) 

103 

39 

3 

141 

15 

Rate  (Cents  Per  Ton) 

357 

63 

239 

330 

133 

Handling  Cost 
(Cents  Per  Ton) 

40 

19 

22 

32 

17 

Rate  (Mills  Per  Ton  Mile) 

228 

53 

464 

107 

198 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 


The  discriminant  analysis  brings  each  variable  into 
the  equation  in  a  stepwise  procedure  based  on  the  contribu¬ 
tion  that  each  variable  makes  to  the  discrimination  between 
groups.  In  this  case  the  order  of  entry  was: 

Step  F  value  to  enter  or  remove 


1.  Rate  137.1935 

2.  Time  of  transit  109.5812 

3.  Length  of  haul  46.3572 

4.  Shipment  size  22.3648 

5.  Size  of  annual  shipment  5.9386 

6.  Handling  costs  1.7916 


The  discriminant  functions  for  each  mode  are  summarized 


below: 

Mcce 

No.  of 
Shipments 

Constant 

Per  Year 

Miles 

Hours 

Shipment 

Sice 

Rate 

Har.vi’  ir.g 
Cost 

Rail 

257 

-  106. 6926a 

3.85272 

-  7.25016 

-  .65165 

4.32353 

31 . 2C829 

S . "So- 7 

3arge 

68 

-  -79.73727 

3.72078 

-  1.61190 

-  .93323 

-.93595 

20.12737 

8.52578 

Truck 

25 

-  -97.00450 

3.22916 

-  5.03098 

-5.93570 

3. 0-5-4 

3 . 2**  5h4 

7.74773 

Rail  Barge 

23 

-  ill. 65173 

4.22C60 

-  5.04177 

-  .68821 

-.31-22 

29.1C273 

3.97553 

Unit  Train 

18 

-  101. 78012 

4.65906 

-  5.72404 

-2.14519 

4.54454 

28,90275 

8.68715 
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In  the  classification  stage  of  the  discriminant  analysis 
each  observation  has  a  value  computed  by  each  of  the  func¬ 
tions.  The  observation  is  assigned  to  the  group  for  which 
the  computed  value  is  nearest  to  zero.  For  example,  com¬ 
puted  discriminant  function  values  for  two  actual  barge 
cases  are  shown  below.  The  first  was  assigned  to  the 
rail  group,  the  second  was  assigned  to  the  barge  group. 


Rail 

Barge 

Truck 

Rail  Barge  Unit  Train 

Case 

1 

1.958 

26.817 

20.839 

7.96 

15.601 

Case 

2 

27.058 

3.782 

57.706 

14.926 

17.301 

Actual 

data  for  each  case 

were : 

Tons 

Per  Year 

Miles 

Hours 

Shipment 

Size 

Rate  Handling 

Case 

1 

2,530 

250 

204 

150 

$5.48  $0,550 

Case 

2 

1,754,000 

309 

108 

11,500 

$  .88  $  .088 

The  data  for  Case  1  looks  like  a  rail  shipment  because  of 
the  shipment  size  and  rate. 

The  next  step  in  estimating  the  demand  function 
is  to  increment  the  rate  variable  for  one  mode  and  reclassi 
fy  the  cases  for  that  mode  by  the  discriminant  functions 
derived  above.  Barge  rates  were  incremented  from  an 
average  of  $  .63  per  ton  to  $1.63  per  ton.  Each  barge 
shipment  was  adjusted  upward  by  $1.00  per  ton.  Figure 
5.1  shows  the  demand  function  derived  from  that  calculation 
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ambiguity,  since  shipment  size  is  a  function  of  which  mode 
completes  the  shipment. 

In  summary,  the  demand  function  estimated  from  all 
coal  shipments  (N=395)  is  price  inelastic  at  current  rates 
and  quantities.  Elasticity  of  price,  estimated  at  current 
rates  and  quantities  is  .62,  that  is  a  one  percent  increase 
in  rate  results  in  a  .62  percent  decrease  in  quantity 
shipped  by  barge. 

Destination  Level.  Observations  were  coded  to  indi¬ 
cate  the  Bureau  of  Economic  Analysis  Region  (BEAR)  of 
destination.  Distinct  characteristics  of  the  demand  func¬ 
tion  for  barge  can  be  noted  at  this  level  of  aggregation. 

Pittsburg  Region.  Table  5.3  shows  the  data  for  35 
coal  shipments  destined  for  the  Pittsburg  region.  Total 
shipments  in  excess  of  44  million  tons  are  hauled  by  three 
modes,  although  barge  accounts  for  98  percent  of  the  traffic. 

TABLE  5.3 

Summary  of  Coal  Shipments  to  Pittsburg  Region 


Rail 

Barge 

Rail/Barge 

No.  of  Shipments 

6 

27 

2 

Annual  Shipment  (1000  Tons) 

39 

1,621 

246 

Average  Shipment  (Tons) 

215 

4,696 

3,950 

Length  of  Haul  (Miles) 

158 

93 

272 

Time  of  Transit  (Hours) 

35 

50 

102 

Rate  (Cents  Per  Ton) 

444 

53 

322 

Handling  Cost  (Cents  Per  Ton) 

27 

20 

11 

Rate  (Mills  Per  Ton  Mile) 

27 

6 

12 

Source :  U . S .  Army  Engineer  Survey ,  1970  and  1971. 
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The  demand  function,  generated  from  data  summarized 
in  Table  5.3  is  essentially  inelastic  near  the  current 
rate,  with  increasing  price  elasticity  thereafter.  Esti¬ 
mated  price  elasticity  of  demand  for  barge  transport  is 
.0002  at  average  rates  and  quantities  represented  in  the 
sample . 

Huntington-Ashland  Region.  Data  for  26  coal  ship¬ 
ments  destined  for  BEAR  52  are  shown  in  Table  5.4.  Total 
annual  shipments  are  in  excess  of  9  million  tons. 

The  demand  function  for  water  transportation  for 
coal  shipments  generated  from  this  data  and  shown  in  Figure 
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TABLE  5.4 


Summary  of  Coal  Shipments  to  Huntington-Ashland  Region 


No.  of  Shipments 

Annual  Shipment  (1000  Tons) 

Average  Shipment  (Tons) 

Length  of  Haul  (Miles) 

Time  of  Transit  (Hours) 

Rate  (Cents  Per  Ton) 

Handling  Cost  (Cents  Per  Ton) 
Rate  (Mills  Per  Ton  Mile) 


Rail 

Barge 

Rail/Bargi 

15 

3 

8 

247 

1,314 

171 

864 

4,600 

1,675 

131 

170 

145 

91 

38 

92 

318 

90 

244 

52 

30 

41 

24 

5 

17 

An  obvious  reason  for  barge  transport  demand  being  inelastic 
in  this  region  is  that  there  are  many  coal  fields  which 
are  as  accessible  by  water  as  by  rail,  and  the  data  summary 
shows  barge  rates  to  be  20  percent  of  rail  rates  and  there 
is  a  time  advantage  by  barge.  Since  a  stepwise  computa¬ 
tional  model  was  used,  there  is  some  evidence  about  the 
relative  importance  of  each  variable.  The  order  of  entry 
of  the  variables  was : 


Step  entered  F  value 

to  enter  or  remove 

1. 

Rate 

6.6814 

2 . 

Shipment  Size 

4.5075 

3. 

Size  of  Annual  Shipment 

2.0051 

4. 

Length  of  Haul 

1.3795 

5. 

Time  of  Transit 

.  3272 

6. 

Handling  Cost 

.  1437 

Cincinnati  Region.  Three  modes  haul  69  coal  ship¬ 
ments  destined  for  the  Cincinnati  Region  (BEAR  62) .  The 
barge  mode  hauls  77  percent  of  the  total  annual  shipment 
(in  excess  of  11  million  tons).  The  share  of  the  total 
hauled  by  each  mode  is  summarized  in  Table  5.5 
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TABLE  5.5 


Summary  of  Coal  Shipments  to  Cincinnati  Region 


Rail 

Water 

Truck 

No.  of  Shipments 

47 

16 

6 

Annual  Shipment  (1000  Tons) 

53 

547 

18 

Average  Shipment  (Tons) 

537 

6,988 

48 

Length  of  Haul  (Miles) 

200 

302 

66 

Time  of  Transit  (Hours) 

156 

92 

18 

Rate  (Cents  Per  Ton) 

449 

98 

232 

Handling  Cost  (Cents  Per  Ton) 

42 

52 

43 

Rate  (Mills  Per  Ton  Mile) 

22 

3 

35 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 


The  demand  function  generated  from  these  data  is 
shown  in  Figure  5.4  Note  that  at  the  current  rate  the 
demand  function  is  price  inelastic,  but  that  the  upper 
range  of  rate  increases  would  result  in  a  price  elastic 
demand  function. 


Louisville  Region.  Data  on  68  shipments  destined 


to  Louisville  Region  are  summarized  in  Table  5.6. 

TABLE  5.6 

Summary  of  Coal  Shipments  to  Louisville  Region 


Rail 

Water 

Truck 

Unit  Train 

No.  of  Shipments 

40 

16 

9 

3 

Annual  Shipment  (1000  Tons) 

26 

290 

3 

984 

Average  Shipment  (Tons) 

744 

5,432 

37 

4,167 

Length  of  Haul  (Miles) 

457 

703 

45 

55 

Time  of  Transit  (Hours) 

438 

168 

5 

55 

Rate  (Cents  Per  Ton) 

715 

183 

3,133 

143 

Handling  Cost  (Cents  Per  Ton) 

38 

77 

109 

23 

Rate  (Mills  Per  Ton  Mile) 

16 

3 

85 

14 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 


The  demand  function  for  water  transportation  gener¬ 
ated  from  these  data  is  shown  in  Figure  5.5.  The  demand 
function  is  essentially  price  inelastic  near  current  price- 
quantity  relationships  but  is  price  elastic  at  higher 
barge  rates.  A  uemand  function  for  rail  transportation 
was  also  estimated  from  the  model.  It  also  shows  price 
inelastic  behavior  near  current  prices,  which  seems  to 
indicate  a  stable  equilibrium  in  intermodal  competition. 

One  important  factor  is  that  rail  hauls  (conventional 

and  unit  train)  are  far  shorter  than  water  hauls.  Shorter 
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Electrical  Utility  Industry.  One  hundred  twenty- two 


Figure  5.6  shows  the  estimated  demand  function  for 
barge  transportation  by  utilities.  The  calculation  under¬ 
estimated  actual  shipments  at  current  prices  by  about 
14  million  tons. 

Metallurgical  Coal.  Coal  shipped  to  iron  and  steel¬ 
makers  is  frequently  mined  in  captive  mines  and  shipped 
on  captive  barges  or  rail  cars.  Captive  mines  and  barges 
are  owned  by  the  coal  and  steelmakers.  The  data  is  sum¬ 
marized  in  Table  5.7. 

TABLE  5.7 

Summary  of  Metallurgical  Coal  Shipments 


Rail 

Barge 

Rail/Barge 

No.  of  Shipments 

31 

9 

17 

Annual  Shipments  (1000  Tons) 

171 

596 

288 

Average  Shipment  (Tons) 

1,183 

4,989 

1,976 

Length  of  Haul  (Miles) 

201 

157 

775 

Time  of  Transit  (Hours) 

107 

113 

233 

Rate  (Cents  Per  Ton) 

419 

74 

398 

Handling  Cost  (Cents  Per  Ton) 

28 

34 

38 

Rate  (Mills  Per  Ton) 

21 

5 

5 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 

The  demand  function,  shown  in  Figure  5.7,  for  barge 
transportation  of  coal  delivered  to  the  iron  and  steel 
industry,  is  somewhat  less  price  elastic  than  of  the  elec¬ 
trical  utility  industry.  Figure  5.7  shows  the  savings 


RATE  C/TON 


resulting  from  barge  shipment  over  shipment  by  rail.  In 
this  case,  benefits  computed  by  traditional  methods,  i.e., 
savings  in  cost  to  shippers,  would  substantially  exceed 
those  computed  by  measuring  the  area  under  the  demand  func¬ 
tion,  since  the  area  ABDE,  representing  savings  computed 
on  average  savings  to  barge  users,  is  greater  than  the  area 
AFE ,  representing  the  area  under  the  demand  schedule  for 
barge  transportation  (or  "willingness  to  pay") . 

Commodities  Other  Than  Coal.  For  all  commodities 
other  than  coal,  the  demand  for  barge  transportation  in¬ 
dicates  the  demand  is  essentially  price  inelastic  at  the 
aggregate  level.  Figure  5.8  shows  the  estimated  demand 
function  based  on  the  total  of  212  shipments  which  include 
data  from  the  Arkansas  River  area.  Table  5.8  shows  the 
data  for  all  commodities  other  than  coal. 
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TABLE  5.8 


Summary  of  Shipments  Other  Than  Coal 


Rail 

Barge 

Truck 

Pipeline 

Pipeline 

Barge 

No.  of  Shipments 

80 

81 

44 

4 

3 

Annual  Shipments  (Tons)  23 

155 

5 

225 

1,973 

Average  Shipment 
(Tons ) 

236 

2,930 

26 

2,354 

12,000 

Length  of  Haul 
(Miles ) 

562 

952 

63 

162 

1,048 

Time  of  Transit 
(Hours ) 

161 

198 

7 

109 

536 

Rate  (Cents  Per  Ton) 

2,415 

283 

901 

120 

202 

Handling  Cost 
(Cents  Per  Ton) 

43 

41 

50 

8 

17 

Rate  (Mills  Per  Ton) 

43 

3 

14 

7 

3 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971 

Chemicals .  For  shipment  of  chemicals,  the  demand 
function  for  barge  transportation  (see  Figure  5.9)  is 
price  inelastic  through  a  substantial  range,  then  highly 
elastic.  One  of  the  established  characteristics  of  chemi¬ 
cal  transportation  is  the  need  for  highly  specialized, 
often  dedicated,  equipment  because  of  corrosiveness  and 
other  problems.  Summary  data  are  presented  in  Table  5.9. 
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TABLE  5.9 


Summary  of  Chemical  Shipments 


Rail 

Barge 

No.  of  Shipments 

3 

10 

Annual  Shipments  (1000  Tons) 

2 

116 

Average  Shipment  (Tons) 

78 

2,924 

Length  of  Haul  (Miles) 

1,100 

981 

Time  of  Transit  (Hours) 

164 

263 

Rate  (Cents  Per  Ton) 

1,733 

296 

Handling  Costs  (Cents  Per  Ton) 

73 

6 

Rate  (Mills  Per  Ton) 

16 

3 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971 


Refined  Petroleum  Products.  Data  on  172  shipments 
were  obtained  in  the  survey  and  are  shown  in  Table  5.10. 

TABLE  5.10 

Summary  of  Shipments  of  Refined  Petroleum 


Rail 

Barge 

Truck 

No.  of  Shipments 

71 

57 

44 

Annual  Shipments  (1000  Tons) 

22 

116 

5 

Average  Shipment  (Tons) 

75 

3,254 

26 

Length  of  Haul  (Miles) 

550 

1,015 

63 

Time  of  Transit  (Hours) 

165 

195 

7 

Rate  (Cents  Per  Ton) 

1,587 

264 

901 

Handling  Costs  (Cents  Per  Ton) 

42 

51 

50 

Rate  (Mills  Per  Ton) 

47 

3 

14 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 


The  demand  function  for  barge  transportation  carry¬ 
ing  refined  petroleum  products  is  price  inelastic.  The 
explanation  is  that  refined  petroleum  products  are  often 
transported  by  user-owned  barges  and  towboats.  Obviously, 
a  huge  rate  differential  exists  between  rail  or  truck 
modes  and  barge.  Large  shipment  sizes  over  long  distances 
are  typical  for  barge. 

Efficiency  of  the  Discriminant  Function.  Besides 
using  the  discriminant  analysis  to  estimate  a  demand  func¬ 
tion,  a  more  obvious  use  in  transportation  analysis  would 
be  to  classify  shipments  with  certain  attributes  by  the 
mode  most  likely  to  be  chosen  by  those  who  make  mode  choice. 

The  discriminant  analysis  program  produces  a  summary 
of  the  results  of  classification  of  the  cases  used  to 
calibrate  the  discriminant  functions.  Obviously  there 
should  be  an  upward  bias  in  the  number  correctly  classified 
because  the  data  are  used  to  calibrate  the  model. ^ 

To  estimate  the  bias  a  run  of  the  model  on  coal 
to  all  destinations  was  made  with  10  percent  of  the  obser¬ 
vations  held  out.  The  hold  out  members  of  the  population 
were  selected  randomly.  The  following  results  were  ob¬ 
tained  : 


^■Ronald  E.  Frank,  et  al,  "Bias  in  Multiple  Discrimi¬ 
nant  Analysis,"  Journal  of  Marketing  Research,  Vol.  II, 
Aug.  1965,  pp.  250-258. 
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Classification  of  380  Observations 


Rail 

Barge 

Truck 

Rail-Barge 

Unit  Train 

1. 

Rail 

175 

10 

38 

30 

2 

2. 

Barge 

1 

36 

8 

2 

16 

3. 

Truck 

1 

0 

21 

0 

0 

4. 

Rail-Barge 

10 

2 

1 

8 

2 

5. 

Unit  Train 

3 

2 

3 

1 

9 

Over  65  percent  of  the  shipments  were  classified  correctly. 
The  classification  functions  were  utilized  to  classify 
holdout  observations  as  follows: 


Classification  of  40 

Holdout 

Observations 

Rail  Barge 

Truck  Rail 

-Barge 

Unit  Train 

1. 

Rail 

21  1 

5 

0 

0 

2. 

Barge 

0  4 

1 

0 

0 

3. 

Truck 

0  0 

1 

0 

0 

4. 

Rail-Barge 

1  1 

1 

1 

0 

5. 

Unit  Train 

0  0 

2 

0 

1 

Of  the  40 

holdout  observations, 

27 

or  68 

percent 

were  classified  correctly. 

Frank  et  al,  suggest  the  following  test  of  signifi¬ 
cance  between  the  proportion  of  correct  classifications 

2 

and  the  proportion  which  could  be  attributed  to  chance: 

2 Ibid,  p.  253. 
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where  Q  =  the  proportion  of  the  sample  observations  correct 
ly  classified  by  the  discriminant  analysis. 

P  =  the  proportion  one  would  expect  by  chance  (if 
groups  are  equal  size,  P  =  .5) 

a  =JEkU 
P  T  n 

t  =  Students  t  distribution 

The  value  of  t  computed  from  the  sample  of  380  cases  is 
5.69,  assuming  the  proportion  one  would  expect  by  chance 
to  be  .5,  and  the  t  value  for  the  40  case  sample  is  2.28. 
Both  estimates  are  significant  at  the  0.05  level  and  there¬ 
fore  there  is  no  evidence  of  significant  bias  due  to  samp¬ 
ling  errors. 

Table  5.11  presents  a  summary  of  classification 
results  for  the  runs  discussed  in  the  demand  analysis 
section. 


Summary 

Demand  functions  for  barge  transportation  for  coal 


and  other  commodities  at  various  levels  of  disaggregation 
have  been  estimated.  Classification  results  were  in  excess 
of  92  percent  correctly  classified  except  for  shipments  of 


refined  petroleum  products.  However,  comparisons  by  ton¬ 
nage  showed  as  much  as  14  percent  of  coal  shipments  by 
barge  to  electrical  utility  users  were  misclassif ied , 
an  indication  that  large  annual  shipments  may  be  more 
likely  to  be  misclassif ied ,  especially  where  unit  train 
technology  is  competing.  Using  a  holdout  sample  of  40 
shipments,  selected  at  random,  the  expected  upward  bias 
in  classification  did  not  materialize.  The  following 
chapter  develops  a  summary  of,  and  conclusions  reached, 
in  this  dissertation. 


TABLE  5.11 


Summary  of  Classification  Results 


%  Correctly 

Demand  Model  Classified 


All  Coal 

95 

Coal  to  Pittsburg  Region 

Destination  BEAR  66 

94 

Coal  to  Huntington-Ashland  Region 
Destination  BEAR  62 

96 

Coal  to  Cincinnati  Region 

Destination  BEAR  54 

97 

Coal  to  Louisville  Region 

Destination  BEAR  52 

93 

Coal  to  Electrical  Utility  Users 

95 

Coal  to  Iron  and  Steel  Industry 

92 

Commodities  Other  Than  Coal 

93 

Refined  Petroleum  Products 

83 

Chemical  Products 

100 

CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 

Characteristics  of  the  Data 

Data  on  nine  transport  modes,  14  major  commodity- 
groups,  totaling  145.9  million  freight  tons  were  analyzed 
for  this  study.  The  number  of  observations  for  automotive 
equipment,  transportation  equipment,  machinery  (except 
electrical) ,  and  furniture  industries  were  too  few  to 
calibrate  a  mode  choice  model.  In  several  commodity  groups 
too  few  observations  were  obtained  for  one  or  more  modes. 
Commodity  groups  for  which  adequate  numbers  of  observa¬ 
tions  were  available  for  calibrating  the  modal  split  model 
were:  coal,  chemicals,  petroleum  products,  primary  metals 

and  fabricated  metals. 

The  data  for  each  mode  and  commodity  were  analyzed 
for  statistical  properties.  In  general,  the  raw  data 
were  found  to  be  skewed  positively  and  characterized  by 
significant  kurtosis.  These  characteristics  diverge  from 
the  assumption  that  the  groups  are  multivariate  normal. 

The  data  were  transformed  to  logarithmic,  exponential 
and  square  root  forms  and  evaluated.  The  logarithmic 
transform  is  characterized  by  small,  negative  skewness 


and  small  values  of  kurtosis.  Thus  the  data  in  logarith¬ 
mic  transform  approach  properties  of  the  normal  distribu¬ 
tion,  at  least  in  the  univariate  case. 


Application  of  Discriminant  Analysis 

Discriminant  analysis  is  based  on  three  basic  assump¬ 
tions  : 

1.  The  probability  density  function  of  the  linear 
discriminant  is  multivariant  normal. 

2.  The  variance-covariance  matrices  for  the  groups 
are  equal. 

3.  The  prior  probability  of  being  in  any  group 
is  equal. 

Watson  reached  the  conclusion  that  since  discriminant 
models  do  not  estimate  probability  of  choice  of  mode  direct¬ 
ly  and  since  discriminant  models  performed  poorly  in  pre¬ 
diction  of  mode  choice  and  classification  that  the  probit 
model  would  serve  the  purposes  of  investigating  the  value 
of  time  in  passenger  traffic.^-  He  found  linear  regression 
gave  an  unbiased  but  inefficient  estimate  of  probability. 

The  SPSS  and  BMD  versions  of  discriminant  analysis 
allow  an  adjustment  for  size  of  groups  which  corrects 
bias  in  classification  due  to  unequal  size  groups.  It  has 
been  shown  that  the  univariant  distribution  of  each 

^Peter  L.  Watson,  The  Value  of  Time:  Behavioral 
Models  of  Model  Choice  (Lexington  Books,  1974),  Chapter  4. 


parameter  approaches  normality  when  transformed  to  logarith¬ 
mic  form.  Box's  M  and  its  associated  F  test  for  equality 
of  group  covariance  is  available  in  SPSS,  but  requires  a 
substantial  increase  in  required  core  storage  space.  Be¬ 
cause  of  limitations  in  core  storage  space  in  the  CDC 
6600  facility  used  for  this  analysis,  the  test  could  not 
be  run.  The  discriminant  model  was  run  in  many  combina¬ 
tions  of  disaggregation  of  shipments  to  explore  the  robust¬ 
ness  of  the  model  under  different  assumptions.  The  model 
shows  a  high  percent  of  correctly  classified  observations 
(from  83  to  100  percent)  which  improves  as  the  classes 
of  shipments  are  disaggregated  to  more  homogeneous  group¬ 
ings.  Tests  of  the  model  using  holdout  data  showed  no 
substantial  upward  bias  in  classification.  Therefore, 
it  is  concluded  that  the  discriminant  models  replicate  the 
behavior  of  freight  transport  users  in  the  Ohio  and 
Arkansas  River  areas . 

Would  alternative  procedures  offer  preferable  re¬ 
sults?  Modal  split  models  calibrated  by  probit  and  logit 
regression  models  offer  the  advantage  of  relaxing  the 
linear  assumption  for  parameters  and  the  assumption  of 
equal  variance/covariance  matrices. 

Oum  criticizes  the  use  of  linear  logit  models  be¬ 
cause  they  impose  many  rigid  a  priori  restrictions  on 
parameters  of  price  responsiveness  of  demand,  including 


l  is 


o 

the  elasticity  of  substitution  and  cross  price  elasticity. 

He  concludes  that  the  structure  of  technology  (or  prefer¬ 
ence)  underlying  the  linear  logit  model  is  severely 
irregular  and  inconsistent.  In  the  price  ratio  versions 
of  the  linear  logit  model,  elasticity  of  substitution 
depends  upon  the  choice  of  the  base  mode  and  the  multi¬ 
nomial  case  restricts  cross  price  elasticities,  with  re¬ 
spect  to  any  given  nonbase  modes,  to  be  identical.  In 
the  price  ratio  case,  the  elasticity  of  substitution  be¬ 
tween  any  two  modes  at  a  given  data  point  is  a  function 
of  the  selected  mode,  however,  the  underlying  elasticities 
of  substitution  do  not  depend  upon  the  selection  of  the 
base  mode.  The  model  does  hold  all  cross  price  elastici¬ 
ties  between  nonbased  modes  to  be  identical.  Therefore, 
while  the  linear  logit  model  may  improve  the  predictive 
ability  of  a  mode  choice  model,  its  utility  to  evaluate 
cross  mode  elasticities  of  substitution  is  severely  limited. 
The  probit  model  contains  similar  problems.  Computational 
capacity  to  calibrate  the  probit  model  by  maximum  likeli- 
hood  methods  have  been  discussed  in  the  literature. 

2 

Oum,  "A  Warning  on  the  Use  of  Linear  Logit  Models 
in  Transport  Mode  Choice  Studies,  "  p.  374. 

^M.J.D.  Powell,  "An  Efficient  Method  for  Finding  the 
Maximum  of  Function  of  Several  Variables  Without  Calcula¬ 
ting  Derivatives,"  Computer  Journal,  1964,  pp .  155-172; 

S.M.  Golfeld,  R.E.  Quandt ,  and  H.F.  Trotter,  "Maximization 
by  Quadratic  Hill  Climbing,"  Econometrica ,  1966,  pp .  541- 
551. 
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However,  it  was  not  possible  to  utilize  an  operational 
version  of  the  maximum  likelihood  methods  for  this  classi¬ 
fication,  since  SPSS  and  BMD  program  packages  do  not 
include  the  method. 

Economic  Interpretation  of  the  Demand  Model 

Even  though  the  statistical  problems  with  the  dis¬ 
criminant  model  appear  to  be  relatively  minor,  some  of 
the  analysis  results  are  difficult  to  justify  economically. 
Some  of  the  simulations  suggest  that  the  demand  for  barge 
transportation  is  significantly  price  inelastic  in  the 
short  run.  Escalation  of  the  price  variable  from  3  to 
5  times  current  relative  prices  are  required  to  shift 
the  slope  to  a  relative  price  elastic  position.  This 
suggests  less  substitution  between  barge  and  other  modes 
than  one  would  expect.  Following  is  a  summary  of  the 


real  price  and  service  ratios  between  barge  and  other 
modes  for  a  number  of  destination  and  commodity  groups: 


TABLE  6.1 


Comparison  of  Cose  and  Service 

3. 

Ratios  Between  Competing  Modes 


Rail 

Barge 


Ratio 


Truck 

Barge 


Ratio 


Unit  Train 
Barge 


Ratio 


All  coal 

rate 

3.7 

7.6  2.85 

speed 

.54 

4.0  .95 

Coal  Dest 

66 

rate 

2.0 

speed 

.52 

Dest 

52 

rate 

4.17 

speed 

.32 

Dest 

62 

rate 

5.2 

9.1 

speed 

.87 

1.15 

Utilities 

rate 

3.2 

speed 

.74 

Chemicals 

rate 

58.7 

speed 

1.79 

Ref.  Petro 

rate 

15.9 

5.0 

speed 

.64 

1.7 

aCost  =  rate 

and  handling  per 

ton  mile;  speed  =  miles  per 

hour 

If  the 

data  are 

representative,  barges  hold  an  un- 

usual  advantage  in  rates  and  a  frequent  advantage  in  speed 
(a  surrogate  for  service) .  The  following  is  a  summary  of 
ratios  of  annual  shipment  and  average  shipment  sizes  for 
various  commodity  destinations. 


TABLE  6.2 


Comparison  of  Annual  Shipment  Size 
and  Average  Shipment  Size  By  Mode 


Rail 

Barge 

Truck 

Barge 

Unit  Train 
Barge 

All  Coal 

Annual 

.  786 

.27 

1.39 

Average 

.197 

.023 

.58 

Coal  Dest  66 

Annual 

.053 

Average 

.  378 

Coal  Dest  52 

Annual 

.188 

Average 

.188 

Coal  Dest  62 

Annual 

.096 

.033 

Average 

.077 

.  007 

Utilities 

Annual 

.  286 

.020 

1.154 

Average 

.  237 

.055 

.53 

Chemicals 

Annual 

.016 

Average 

.027 

Ref.  Petro 

Annual 

.191 

.04 

Average 

.023 

.007 

In  all 

cases  barge  users 

are 

large  consumers 

and  can 

receive  larg 

e  shipment  sizes. 

The 

combined  impact 

of 

lower  rates, 

good  relative  service 

characteristics 

,  and 

propensity  to  serve  large  customers  suggest  the  limited 
economic  substitutability  of  other  transport  modes  given 


the  same  shipment  characteristics  of  barge,  except  when 
other  modes  can  combine  technology  and  cost  effectiveness 
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as  is  the  case  with  pipelines  and  unit  trains.  If  small 
users  are  partitioned  out,  only  unit  train  technology 
and  pipeline  offers  a  substantial  economic  substitution. 

What  does  this  mean  to  the  policy  issue  of  waterway 
user  charges  and  increasing  congestion  at  certain  points 
in  the  waterway  system?  This  analysis  shows  that  user 
charges  are  unlikely  to  displace  much  traffic  from  barge 
to  rail  in  the  short  run.  Sasaki  showed  little  long  run 
mode  substitution  in  coal  shipments  by  electrical  utility 
users. ^  Technological  adjustment  by  the  electrical  utility 
industry  is  characterized  by  long  lead  times  (10-20  years) 
and  thermal  plants  have  a  normal  useful  life  of  30  to 
35  years.  However,  the  location  of  coal  mine  development 
could  shift  toward  areas  served  primarily  by  rail  if  barge 
rates  rise  and  service  deteriorates  due  to  the  waterway 
congestion.  There  are,  however,  countering  trends.  Rail 
rates  for  coal  are  increasing  faster  than  rail  rates  for 
other  commodities.  Unit  train  operations  are  creating 
accelerated  maintenance  problems  in  track  and  roadbed. 
Utilities  use  their  coal  inventory  policies  as  a  hedge 
against  potential  strikes  in  the  coal  industry  and  generally 
prefer  to  keep  part  of  their  coal  acquisition  scattered 

^Sasaki,  A  Regional  Model  of  the  Future  Demand  for 
Transportation:  The  Case  of  Barge  Transportation. 


1 

I 

f 

among  small  independent  mines  who  normally  use  truck  and 
rail.  Therefore,  utility  users  keep  more  transport  tech¬ 
nology  alternatives  open  to  them  and  maintain  3-4  months 
stockpiles  of  coal  to  deal  with  uncertainty  in  coal  pro-  ! 

i 

duction  and  with  potential  service  problems  with  transport  \ 

i 

modes. 

I 

Other  commodity  users  have  special  conditions  which 
influence  their  choice  of  mode  and  demand  for  barge  trans-  ■ 

port.  Chemical  products  are  often  corrosive  and  to  some 
degree  a  safety  hazard.  Users  frequently  own  specialized 
equipment  to  haul  chemicals  whether  by  rail,  truck  or 

barge.  Substitution  is  limited  by  this  technological  ; 

specialization,  especially  in  the  short-run.  Barge  trans- 

i 

port  has  a  safety  advantage  since  the  adverse  impacts 
of  an  accident  can  be  mitigated  by  the  dilution  effects 
of  water.  This  can  backfire  if  water  supply  intakes  are 
located  near  a  spill  site.  However,  in  general,  hazardous 
and  corrosive  chemicals  are  much  less  likely  to  cause 
health  and  safety  problems  if  shipped  by  water. 

The  divergence  of  mode  options  between  large  and 
small  users  are  significant  in  chemicals,  petroleum  pro- 

i 

ducts  and  to  a  lesser  degree,  users  of  iron  and  steel 
products.  Large  users  can  accommodate  barge  load  quanti¬ 
ties  and  use  enough  annually  to  justify  barge  and  dock 
investments.  Small  users  do  not  have  this  option  unless 

I 

i 
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an  intermediary  can  introduce  a  barge  handling  service. 

Most  waterways  induce  the  development  of  wholesaler/ 
jobbers  who  can  handle  barge  load  quantities  and  pass 
along  part  of  the  economic  advantage  of  water  transporta¬ 
tion  to  smaller  users. 

Weakness  of  this  Analysis.  Although  there  was  no 
attempt  to  construct  a  theoretical  and  empirical  analysis 
of  the  long-run  economic  demand  for  barge  transportation, 
most  significant  investment  and  policy  questions  require 
an  anlysis  of  the  long-run.  If  the  commodity  users  of 
barge  transportation  are  confronting  constant  cost  and 
predictable  parametric  shifts  in  the  demand  function, 
and  if  the  relative  prices  of  barge  and  its  competing 
modes  are  expected  to  be  constant,  a  competent  analysis 
of  the  short-run  can  satisfy  a  large  share  of  the  long- 
run  policy  issues.  If  these  conditions  are  not  anticipated, 
a  substantial  number  of  projection  and  analytical  issues 
are  raised.  The  earlier  discussion  of  the  coal  using 
electrical  utility  suggest  that  a  fair  appraisal  can  be 
made  from  short-run  behavior  models  since  the  technological 
state  of  the  industry  can  change  only  very  slowly.  In 
addition,  the  uncertainties  of  determining  the  economic 
and  environmental  efficacy  of  nuclear  power,  the  environ¬ 
mental  policies  of  handling  effluent  requirements  of  coal 
fired  plants  and  the  risk  involved  from  sole  dependence 
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on  softer  energy  producing  technology  suggest  that  an 
enlightened  policy  is  to  keep  many  options  open.  This 
entails  substantial  economic  costs,  and  the  uncertainties 
involved  in  the  best  available  analytical  capability  cannot 
be  reduced  to  provide  an  obvious  limited  set  of  choices 
with  small  risk. 

The  data  acquired  for  this  analysis  would  have  been 
greatly  enhanced  if  more  data  on  service  time  and  handling 
costs  for  individual  shipments  for  each  origin-destination 
pair  had  been  obtained.  More  observations  would  have 
provided  a  basis  for  determining  the  variance  around  the 
expected  service  time  and  therefore  would  have  allowed 
the  explicit  introduction  of  uncertainty  and  service  time 
into  the  modal  split  model.  Much  of  the  rail  data  was 
obtained  at  a  time  when  the  New  York  Central-Pennsylvania 
system  was  operating  at  its  worst.  Thus  relative  service 
times  obtained  for  rail  may  have  been  biased  due  to  below 
normal  performance. 

Handling  costs,  defined  as  those  costs  entailed 
in  excess  of  the  transport  rate  to  complete  delivery  and 
unloading,  are  difficult  to  estimate.  With  rare  exceptions, 
handling  costs  are  not  identified  in  shipper  or  receiver 
cost  accounting.  Thus  the  persons  interviewed  would  have 
to  estimate  the  annual  manpower  and  equipment  costs  in¬ 
volved  in  handling  and  divide  by  the  tonnage  involved  to 
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approximate  the  handling  cost  estimate.  The  data  from 
electrical  utility  users  are  considered  to  be  better  since 
they  make  a  substantial  investment  in  material  handling 
to  unload  coal  and  normally  have  the  costs  identified 
in  their  capital  accounts.  Yet  the  material  handling 
costs  due  to  managing  three  to  four  month  coal  stockpiles 
are  significant  and  certainly  not  solely  attributable 
to  the  transport  mode  selected. 

Transportation  data  seldom  are  obtained  in  a  way 
to  correlate  it  with  economic  trade  flows .  Although  barge 
shipments  are  monitored  in  a  complete  annual  census  by 
the  Corps  of  Engineers,  rail  shipment  data  are  available 
only  from  a  one  percent  way  bill  sample  collected  by  the 
Interstate  Commerce  Commission  and  systematic  national 
truck  traffic  is  unregulated.  There's  no  good  way  to 
know  whether  the  data  collected  for  this  or  any  other 
study  is  consistent  with  actual  trade  flows.  This  problem 
is  receiving  more  attention  at  the  Department  of  Transporta¬ 
tion  and  better  overall  data  are  becoming  available.  How¬ 
ever,  transportation  data  are  greatest  in  number  about 
physical  rather  than  the  economic  flows. 

The  Need  for  Behavioral  Models  for  Transportation 
Demand.  Transport  planners  tend  to  think  of  the  trans¬ 
portation  as  a  physical  network  flow  problem.  This  is 
because  many  are  trained  in  engineering  or  geography. 
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In  addition,  the  nation's  transport  system  reflects  a 
mixed  public  and  private  decision-making  structure.  Most 
transport  models,  even  economic  ones,  tend  to  be  designed 
under  the  rubric  of  a  social  optimum  which  places  a  single 
monopolistic  decision-maker  in  control."*  As  long  as  indi 
vidual  shippers  control  the  decisions  about  how  much  to 
ship,  by  what  route  and  mode,  there  is  a  substantial  basi 
to  develop  the  demand  models  to  replicate  the  behavior 
of  these  independent  decision-makers.  Aggregate  models 
tend  to  ignore  or  to  place  different  weights  on  perceived 
costs  of  elements  of  transport  demand  than  the  individual 
decision-makers.  Aggregate  models  tend  to  ignore  or  to 
place  different  weights  on  perceived  costs  of  elements 
of  transport  demand  than  the  individual  decisions.  The 
best  example  of  this  tendency  is  the  propensity  to  use 
rates  as  a  basic  and  frequently  only  element  of  a  cost 
minimizing  transport  flow  model. ^  The  consideration  of 
nonpriced  service  elements ,  the  tendency  to  handle  tech¬ 
nology  as  an  independent  factor,  and  the  tendency  to  use 
aggregate  cost  minimizing  models  bias  the  actual  behavior 

■*Howe  et  al .  ,  Inland  Waterway  Transportation: 
Studies  in  Public  and  Private  Management  and  Investment 
Decisions,  see  Chapter  1. 

^CACI,  Inc.,  Potential  Impacts  of  Selected  Inland 
Waterway  User  Charges,  U.S.  Army  Corps  of  Engineers, 
Washington,  D.C.,  December  1976. 
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of  transport  users.  Demands,  or  requests  from  transporta¬ 
tion  users  and  other  interested  participants  in  public 
transport  policy  formulation,  for  a  more  realistic  analy¬ 
sis  can  be  facilitated  by  an  emphasis  on  disaggregate 
behavioral  models.  The  data  needs  are  little  more  than 
would  be  required  to  calibrate  more  highly  aggregated 
models,  and  there  is  a  variety  of  analytical  procedures 
available  to  calibrate  the  models. 
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APPENDIX  A 


DESCRIPTION  OF  DATA 

This  appendix  displays  and  critiques  the  data  used 
in  this  study.  Basic  data  gathering  was  performed  by 
a  U.S.  Army  Engineer  Division,  Ohio  team  during  the  summers 
of  1970  and  1971  and  a  U.S.  Army  Engineer  Division,  South¬ 
western  team  during  1971.  A  total  of  815  observations 
are  available  with  complete  data  for  the  attributes  of 
shipments  and  mode  characteristics  needed  for  the  mode 
choice  model: 

Quantities  shipped  per  year  in  tons 

Miles  of  haul 

Time  of  transit  (hour) 

Shipment  size  (ton) 

Transport  rate  (cents  per  ton) 

Handling  costs  (cents  per  ton) 

The  observations  were  coded  for  commodity  groups 
by  SIC  codes  at  the  4  digit  level.  All  analyses  group 
the  commodities  at  the  2  digit  level  as  shown  in  Table  A.l. 


TABLE  A . 1 


Commodity  Groups  by  SIC 
and  Number  of  Observations 


Commodity  Group 

SIC  Code  No.  of 

Observations 

Coal 

12 

398 

Crude  Petroleum 

13 

16 

Food  (grains) 

20 

22 

Furniture 

25 

2 

Paper 

26 

28 

Chemicals 

28 

27 

Refined  Petroleum  Products 

29 

177 

Rubber  and  Plastic  Products 

30 

12 

Primary  Metal  Products 

33 

48 

Fabricated  Metal  Products 

34 

64 

Machinery,  except  electrical 

35 

5 

Electrical  Machinery 

36 

14 

Transportation  Equipment 

37 

1 

Automotive  Equipment 

50 

2 

Source:  U.S.  Army  Engineer 

Surveys,  1970  and 

1971. 
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Table  A. 2  summarizes  the  data  by  commodity  group, 
average  tons  per  year  and  total  tonnage  per  year  by  mode. 

TABLE  A. 2 

Summary  of  Sample  Data 
Tons  Per  Year  by  Commodity  by  Mode, 

1970-1971 


Commodity  Group 

Number  of 
Observations 

Average  Tons 
Per  Shipment 
Per  Year 

Total 
Tons  Per 
Year  (000) 

Coal 

398 

307,141 

122,203 

Rail 

257 

81,141 

20,853 

Barge 

68 

1,013,449 

68,915 

Truck 

26 

27,418 

713 

Rail-Barge 

23 

262,840 

6,045 

Truck-Barge 

6 

15,449 

93 

Truck-Rail -Barge 

7 

8,000 

56 

Unit  Train 

18 

1,418,222 

25,528 

Crude  Petroleum 

16 

702,712 

11,24/' 

Rail 

1 

72,800 

73 

Barge 

10 

466,098 

4,661 

Pipeline 

2 

295,887 

592 

Pipeline-Barge 

3 

1,972,610 

5,918 

Grain 

22 

6,315 

135 

Rail 

18 

5,379 

97 

Barge 

3 

13,708 

41 

Truck 

1 

1,000 

1 

Furniture 

2 

1,800 

3.6 

Rail 

1 

360 

.4 

Barge 

1 

3,240 

3.2 

Paper 

28 

3,714 

104 

Rail 

8 

7,625 

61 

Truck 

16 

2,687 

43 
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TABLE  A. 2  -  Continued 


Number  of 

Commodity  Group  Observations 

Average  Tons 
Per  Shipment 
Per  Year 

Total 

Tons  Per 
Year  (000) 

Chemicals 

27 

46,822 

1,265 

Rail 

9 

5,742 

52 

Barge 

13 

92,425 

1,202 

Truck 

5 

2,200 

11 

Petroleum 

Products 

177 

49,202 

8,709 

Rail 

73 

21,941 

1,602 

Barge 

57 

115,611 

6,590 

Truck 

44 

4,638 

204 

Rail-Barge 

1 

2,725 

3 

Pipeline 

2 

155,203 

310 

Rubber  and 

Plastic  Products 

12 

1,572 

17 

Rail 

3 

863 

3 

Truck 

9 

1,803 

16 

Primary  Metals 

48 

38,477 

1,846 

Rail 

27 

59,202 

1,598 

Barge 

6 

30,167 

181 

Truck 

15 

4,496 

67 

Fabricated  Metals 

64 

3,745 

240 

Rail 

25 

1,563 

39 

Barge 

7 

22,421 

157 

Truck 

20 

1,727 

35 

Rail-Barge 

1 

2,000 

2 

Truck-Barge 

11 

647 

7 

Machinery 

Except  Electrical 

4 

1,995 

8 

Rail 

? 

1,742 

3 

Barge 

1 

3,996 

4 

Truck-Barge 

1 

500 

1 

l  34 


1 


TABLE  A. 2  -  Continued 


Commodity  Group 

Number  of 
Observations 

Average  Tons 
Per  Shipment 
Per  Year 

Total 

Tons  Per 
Year  (000) 

Electrical  Machinery  14 

8,851 

126 

Rail 

1 

15,000 

15 

Barge 

6 

11,833 

71 

Truck 

5 

782 

6 

Truck-Barge 

2 

17,000 

34 

Transportation 

Equipment 

1 

5,000 

5 

Rail 

1 

5,000 

5 

Wholesale  Trade 
(Automotive 

Equipment) 

2 

28,500 

57 

Rail 

2 

28,500 

57 

TOTAL 

815 

145,963 

Source:  U.S.  Army  Engineer  Survey,  1970  and  1971. 


For  595  observations  there  is  complete  origin  and 
destination  coding.  Table  A. 3  shows  the  number  of  observa¬ 
tions  by  origin  and  destination  for  the  data  from  the 
Ohio  and  Arkansas  River  areas.  The  origin  and  destination 
code  is  the  designation  of  the  Department  of  Commerce 
for  functional  economic  areas  (see  Figure  A.l),  Table  A. 4 
summarizes  the  total  data  set  by  mode  and  attributes  of 
modal  choice. 
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Table  A. 5  shows  cost  per  ton  mile  and  speed  for  the 
entire  sample  by  mode.  Speed  was  computed  by  dividing 
miles  of  transit  by  hours  in  transit. 

Since  the  form  of  the  modal  split  and  demand  model 
is  sensitive  to  the  distributional  characteristics  of  the 
variables,  these  attributes  were  evaluated  for  each 
variable.  In  particular,  skewness  and  kurtosis  show  how 
well  the  data  fit  the  properties  of  a  normal  distribution. 

In  Table  A. 6  and  Figure  A. 2  the  natural  data  shows  substan¬ 
tial  variation  from  desired  properties  especially  for  annual 
tonnage  and  rate.  The  values  for  data  transformed  into 
natural  logarithms  and  into  the  exponential  form  are  also 
shown.  The  exponential  transform  uniformly  exaggerates 
the  variations  whereas  the  logarithmic  transform  provides 
almost  perfect  properties  of  a  normal  distribution.  The 
logarithmic  transform  performs  much  better  in  discriminant 
and  regression  analyses. 
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COMPARISON  OF  SKEWNESS  i  XURTOSIS  VALUES 
Of  TRANSFORMED  VARIABLES 
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TABLE  A. 5 


Cost  Per  Ton  Mile,  Speed  and 
Length  of  Haul,  All  Commodities 


Mode 

N 

Cost  Per 
Ton  Mile 
(Cents) 

Speed 

(mph) 

1 

Length  of 

Haul  (Miles) 

Rail 

427 

3.07 

2.47 

343 

Barge 

171 

.50 

3.98 

673  ! 

Truck 

145 

5.31 

9.63 

248  u 

t 

Rail-Barge 

25 

.72 

2.84 

617 

, 

Truck- Barge 

16 

1.16 

2.49 

783  f 

Rail-Truck-Barge 

4 

3.57 

.78 

ib»  ;1 

Unit  Train 

18 

1.71 

2.56 

117 

Pipeline 

4 

.79 

1.49 

162 

Pipeline-Barge 

3 

.21 

1.95 

1,048  j 

Summary  of  Data  by  Variable 


